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The  objective  of  the  effort  was  to  experimentally  determine  the 
influence  of  the  atmosphere  on  radar  resolution  in  range,  angle  and  velocity. 
The  approach  was  to  utilize  satellite  borne  transmitters  and  ground  based 
receivers  designed  such  that  measurements  could  be  made  of  (1)  the  spectrum 
of  phase  fluctuations,  (2)  the  angular  spectrum,  and  (3)  dispersion. 


Technical  and  financial  problems  negated  the  possibility  of  conducting 
the  satellite  experiment.  Radio  astronomy  experiments  were  performed  howeve 
to  deduce  the  phase  fluctuations  contributed  by  the  ionosphere  and  aircraft 
borne  sources  were  utilized  to  deduce  the  phase  fluctuations  induced  by  the 
troposphere.  The  tropospheric  measurement  of  phase  fluctuations  was  limited 
by  the  inability  to  accurately  determine  aircraft  position;  therefore  the 
resultant  spectrum  of  phase  fluctuations  was  contaminated  by  aircraft  induced 
phase  variations.  It  is  still  useful  within  limits  however. 


Estimates  of  achievable  resolution  were  made  based  on  data  collected 
by  the  contractor  and  the  National  Bureau  of  Standards.  The  following 
conclusions  were  reached: 


(1)  Angular  resolution  of  0. 1  milliradian  can  be  achieved  at  frequencies 
greater  than  100  mcs  for  targets  of  equal  cross  section.  At  X-band  the  cross 
section  may  vary  as  much  as  13  db  and  still  be  resolved. 

(2)  Velocity  resolution  of  0. 1  ft  per  second  can  be  attained  for  targets 
of  equal  cross  section.  At  X-band,  this  degree  of  resolution  can  be  achieved 
on  targets  whose  cross  section  differs  by  13  db. 

(3)  Range  resolution  of  1,0  foot  can  be  achieved  at  frequencies  greater 
than  4,000  mcs  for  targets  of  equal  cross  section.  Although  difficult  to  evaluate 
it  is  expected  that  a  difference  of  target  cross  section  of  30  db  will  still  permit 
resolution  to  1.0  foot. 

The  above  estimates  are  considered  to  represent  worst  case  conditions; 
therefore  system  designers  can  utilize  these  values  with  confidence.  It  is 
evident  from  this  effort  that  worst  case  estimates  of  the  atmospheric  influence 
on  radar  resolution  can  be  provided  utilizing  indirect  measurement  techniques 
but  that  the  validity  of  these  estimates  become  highly  questionable  when  it  is 
necessary  to  establish  limits.  It  is  conceivable  that,  as  radar  resolution 
requirements  become  more  stringent,  processing  techniques  will  be  required 
which  will  utilize  statistical  descriptions  of  a  signal  in  time  and  space  to  negate 
the  degrading  influence  of  the  medium.  It  has  been  realized  from  this  effort 
that  costly,  sophisticated  approaches  are  necessary  to  obtain  the  required 
statistical  representation. 
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ABSTRACT 


The  results  of  a  theoretical  and  experimental  program  to  determine 
the  influence  of  the  atmosphere  on  radar  resolution  are  presented.  The  limits 
imposed  oy  the  atmosphere  on  resolution  in  range,  angle  and  velocity  are 
treated  theoretically  and  estimated  values  are  presented  based  on  current 
experimental  data.  General  expressions  are  derived  which  permit  assessment 
of  the  atmosr.heric  influence  on  resolution  as  a  function  of  frequency.  Consider¬ 
ation  is  given  to  multipath  dispersion,  phase  fluctuations  and  unequal  target 
cross  sections. 

A  description  of  a  radio  astronomy  experiment  which  was  performed  to 
deduce  the  phase  fluctuations  incurred  by  a  wave  traversing  the  ionosphere  is 
presented  with  the  results  obtained.  Similarly  an  experiment  performed  to 
determine  the  tropospheric  contribution  to  phase  path  fluctuations  is  presented 
along  with  a  discussion  of  the  validity  of  the  data  as  limited  by  the  measurement 
technique  employed.  An  application  of  the  data  is  included  which  provides 
an  assessment  of  the  influence  of  phase  fluctuations  on  the  performance  of  a 
focussed,  multiple  antenna,  communications  system. 
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INTRODUCTION 


The  Wide  Band  Program  was  instituted  to  determine  the  limitations 
imposed  upon  radar  resolution  in  angle,  range  and  velocity  as  a 
consequence  of  propagation  through  the  atmosphere  of  the  earth.  Insofar 
as  the  effects  of  a  nuclear  detonation  upon  radar  resolution  are  represented 
by  naturally  occurring  auroral  phenomena,  the  Wide  Band  program  would 
also  determine  the  effects  of  a  high  altitude  nuclear  detonation  upon 
radar  resolution. 

The  Wide  Band  program  embraces  both  theoretical  and 
experimental  investigations.  The  experimental  investigations  were 
to  include  but  not  necessarily  be  limited  to  the  following  techniques: 

1)  A  passive  satellite  experiment  utilizing  ground>based  trans¬ 
mitters  and  receivers; 

2)  An  active  satellite  experiment  with  ground  based  receivers 
and  satellite -borne  transmitters; 

3)  An  aircraft  experiment  using  airborne  transmitters  and 
ground-based  receivers; 

4)  A  radio  astronomy  experiment  using  an  extra-terrestrial 
radio  source. 
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This  is  the  summary  report  for  the  Wide  Band  Program.  We 
shall  briefly  review  the  theoretical  aspects  of  resolution,  indicate 
those  experimental  techniques  that  were  deemed  desirable  and  feasible 
and  discuss  the  reasons  for  rejecting  certain  of  the  experiments.  The 
above  discussions  will  be  somewhat  brief  since  many  of  these  topics 
have  been  treated  in  some  detail  in  referenced  Wide  Band  reports. 

The  report  will  conclude  with  a  summary  of  our  best  estimates  of  the 
atmospheric  limitations  to  radar  resolution. 


2 


i.  RADAR  RESOLUTION 


Although  th«  general  concepts  of  radar  resolution  are  readily 
comprehended,  a  discussion  of  the  resolution  of  specific  radar  systems 
is  not  at  all  straightforward.  The  general  concepts  of  resolution, 
derived  from  the  classical  optical  "Rayleigh'*  resolution,  criteria 
(Jenkins  and  White.  1950)  define  resolution  as  the  ability  to  distinguish 
between  two  radar  targets  (of  equal  cross  section)  which  difier  in  only 
one  coordinate  angle  (azimuth  or  elevation),  range  or  velocity.  In 
the  cases  of  interest  however,  the  two  targets  may  have  far  difierent 
radar  cross  sections  and  the  radar  echctes  from  each  target  may  be 
time  varying  because  of  target  motion  or  atmospheric  effects.  The 
time  varying  nature  of  the  signal  complicates  the  problem  of  system 
evaluation. 

Before  discussing  the  detailed  relationship  among  resolution, 
radar  system  parameters  and  atmospheric  effects,  it  might  be  useful' 
to  reiterate  the  '  /ell  known  definition  of  radar  cross  section: 


’T  •  Pr 


Pr^rfl  Pr^rP 


■rr  ^  / 


lp,(€)r  V 


radar  cross  section 
received  power 
range  to  target 
transmitted  power 
antenna  gain 

effective  aperture  of  the  antenna 
temporal  variation  of  the  received  signal  voltage 

It  is  important  to  note  that  since  the  cross  section  is  defined 
in  terms  of  tl^e  received  power  (and  the  system  parameters)  there  is 
an  integration  time  required  to  ascertain  the  received  power.  Whenever 
there  is  relative  motion  between  the  target  and  the  radar  system  (due 
to  either  translation  or  rotation  or  both)  the  received  signal  voltage 
will  vary  for  all  targets  except  a  single  point  target.  When  one  takes 
signal-to-noise  ratios  or  propagation  fluctuations  into  account,  the 
received  signal  voltage  for  even  a  single  point  target  will  also  vary 
with  time.  With  no  a  priori  knowledge  of  the  characteristics  of  the 
target,  we  are  required  to  average  the  received  signal  to  determine 
the  cross  section. 

The  preceding  remarks  have  been  included  to  indicate  that  when 
discussing  resolution  in  terms  of  relative  target  "cross  section"  and 
real  radar  systems,  a  certain  amount  of  time  averaging  will  be  required 
for  each  and  every  system.  If  the  averaging  time  required  to  account 


where  - 

Pr  = 

R  = 

A  = 


4 


for  the  relative  motions  of  the  targets  and  the  actual  signaUto-noise 
ratio  is  very  short  compared  to  the  time  necessary  for  the  atmospheric 
effects  to  change  noticeably,  then  the  effects  of  the  atmosphere  upon 
radar  resolution^  may  be  negligible. 

It  would  seem  however,  that  the  integration  time  of  most  systems  is 
such  that  one  should  consider  the  effcxts  of  the  atmosphere  upon  radar 
system  resolution  limitations. 
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2.  ANGULAR  RESOLUTION 


Angular  resolution  is  usually  defined  as  the  ability  of  a  system 
to  resolve  two  targets,  of  equal  cross  section  at  the  same  -ange, 
differing  only  in  their  angular  coordinates.  Note  that  this  definition 
says  nothing  about  the  accuracy  of  location  of  a  single  target.  The 
angular  resolution  of  a  monostatic  radar  system  is  approximately 
equal  to  the  antenna  beamwidth. 

The  antenna  beamwidth  and  its  sidelobe  levels  are  defined  in 
terms  of  the  output  power  of  the  antenna  for  different  orientations 
of  the  antenna  in  the  presence  of  a  simple  plane  wave.  In  a  real 
atmosphere,  due  to  ever-present  -albiet  normally  weak  forward  scattering 
by  inhomogeneities  of  refractive  index, multipath  propagation  can  exist 
and  the  incoming  signals  cannot  be  completely  represented  by  a  single 
pU  ne  wave  but  rather  must  be  represented  by  a  summation  of  plane 
waves  in  angle  space.  Due  to  the  random  nature  of  the  forward  scattering 
process  giving  rise  to  the  multipath  propagation,  the  phase  relationship 
between  the  plane  wave  components  in  the  voltage  angular  spectrum  is 
random  and  it  is  usual  to  speak  of  the  time  averaged  angular  power 
spectrum  of  plane  waves  which  describes  the  power  associated  with  each 
component  of  the  angular  spectrum.  Under  these  circumstances  the 
antenna  is  said  to  act  as  a  space  filter  described  by  the  antenna  power 
pattern  and  the  average  power  output  is  the  integral, over  all  angle  space, 
of  the  product  of  the  antenna  power  pattern  and  the  angular  power  spectrum. 
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The  angular  pouer  spectrum  describes  on  the  average,  the 
distribution  of  power  in  the  received  signal  as  a  function  of  the  ai^[le 
of  arrival.  The  emphasis  given  to  the  term  on  the  average  is  warranted. 
The  angular  power  spectrum  does  not  describe  the  time  history  of 
the  angle  of  arrival  of  radio  waves.  When  assessing  the  effects  of 
diffractive  scattering  upon  resolution,  it  is  important  to  consider 
the  statistics  of  the  incoming  signal  on  the  same  time  basis  as  the 
time  required  to  move  the  antenna  through  an  angular  distance  equal 
to  at  least  two  resolution  elements.  If  one  could  move  the  antenna  through 
this  angular  distance  in  a  time  very  short  compared  to  the  time  constant 
of  the  atmospheric  processes  (say  for  example  10  milliseconds),  one 
might  observe  the  energy  coming  from  a  direction  slightly  different 
from  the  line  of  sight  path  between  the  observer  and  the  target,  but 
the  beam  pattern  of  the  antenna  would  be  preserved.  Two  targets 
separated  by  more  than  the  angular  beamwidth  of  the  antenna  might  be 
resolved  on  the  basis  of  this  short  sample  if  the  target  field  strengths 
were  well  above  noise.  Some  time  later,  the  apparent  angles  of  arrival 
of  the  two  targets  might  be  different  from  before  but  again  the  beam 
pattern  of  the  antenna  might  be  preserved.  Under  these  circumstances 
the  effect  of  multipath  propagation  would  cause  the  instantnaeous  signal 
to  "glint"  over  a  range  of  directions  centered  on  the  line  of  sight  path. 
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So  far  we  have  described  the  principles  which  might  be  used 
to  determine  the  effects  of  multipath  or  forward  scatter  upon  angular 
resoluti^.  We  must  now  qualitatively  assess  these  effects  upon  radar 
resolution. 


2. 1  Angular  Spectrum  and  Angular  Resolution 

Bramley  (1951)  has  derived  the  angular  spectrum  of  plane  waves 
produced  by  a  pliase  changing  screen  in  which  the  spatial  autocorrelation 
function  of  the  fluctuations  of  phase  was  of  the  form 

4  « '  y  ■  (tJ” 

and  found  that  the  angular  power  spectrum  is  given  by  (equation  36  of 
Bramley) 

■ 

where  is  the  mean  squared  phase  fluctuation 

is  the  "correlation  distance"  of  the  phase  fluctuations  measured 
in  wavelengths 

is  the  sine  of  the  angle  of  arrival  measured  from  the  vertical 
is  a  delta  function 
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The  delta  function  in  equation  (2)  repreaenta  the  aero  order 
undeviated  plane  wave  -  that  wave  which  would  propagate  in  the  abaence 
of  the  phaae  fluctuationa.  Aa  the  mean  aquare  phaae  fluctuaticma  increaae, 
power  ia  abatracted  from  the  aero  order  undeviated  wave  and  ia  fed  into 
the  continuoua  apectrum  of  plane  wavea  repreaented  by  the  aecond  term 
of  equation  (2).  One  would  expect;  theiv  that  aa  the  mean  aquared  refractive 
index  fluctuationa  increaae  the  power  in  the  delta  function  in  the 
spectrum  will  decrease  to  a  value  such  that  it  ia  indistinguishable  from 
the  continuous  spectrum  and  the  incoming  spectrum  is  then  spread 
over  a  range  of  angles.  In  Appendix  A  we  show  that  when  ^  is  of  the 
order  of  1. 5  radians  the  power  associated  with  the  delta  function  is 
1/10  of  the  power  associated  with  the  continuoua  spectrum  and  the  angular 
spectrum  has  a  half*width«  /X  d  given  by 

-  .  {,<)  Vo 

><o 

and  that  therefore  the  minimum  angular  resolution  which  can  be  expected 
is  ^  ^ 

(}  2.2  Resolution  Limitations  Imposed  By  The  Auroral  Ionosphere 

An  experiment  to  determine  the  ionospheric  limitations  to  angular 
resolution  was  conducted  under  the  Wide  Band  Program.  This  experiment 
utilized  emissions  from  the  radio  source  in  Cassiopeia  A  and  observed 
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the  phase  fluctuations  intmduced  by  the  auroral  ionosphere  on  frequencies 
of  5>Z  and  147  Mes  simultaneously  on  baselines  of  300,  oOO,  and  <>00  feet 
East-West  and  1000  and  200  feet  North-South.  A  complete  description 
of  the  experimental  technique  and  the  detailed  results  can  be  found  in 
the  Cornell  Aeronautical  Laboratory  Report  No.  UB-1363-P-4,  January 
l<>o2.  A  supplementary  report,  covering  an  additional  six  months  of 
work  is  appended  to  this  report  as  Appendix  B.  To  sum  up  the  results 
of  the  experimental  investigations  of  the  ionosphere,  it  appears  that 


a  radio  wave,  traversing  the  entire  ionosphere  will  suffer  phase 
fluctuations  whose  mean  square  value  Vi  is  given  by  V)  a.  A 


square  radians  and  ,  the  correlation  <listance  of  the  phase  fluctuations 


IS  of  the  order  of  200-300  meters. 


For  frequencies  well  in  excess  of  300  Mes  (A  =  1  rneter)  the 
mean  squared  phase  fluctuations  in  the  ionosphere  will  be  much  less  than  one 
square  radian  and  the  angular  spectrum  of  the  wave  as  received  at 
the  ground  will  consist  of  a  well-defined  angle  of  arrival  with  some  weak 
continuous  scatter.  For  such  frequencies,  the  radar  system  angular 
resolution  will  be  limited  only  by  tne  physical  size  of  the  aperture, 
"Resolution"  in  the  preceding  sentence  is  used  in  the  context  of 
distinguishing  two  targets  of  equal  or  nearly  equal  cross  section. 
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For  frequencies  lees  than  300  Mce  (  ^  i  meterl  ionospheric 

phase  fluctuations  will  set  a  limit  to  the  maximum  obtainable  resolution. 

In  Radio  Astronomy  Report  No.  UB*i363>P>4  we  developed  the  formula 
(derived  a  little  more  elegantly  in  Appendix  A)  ^  where  J  is 

the  correlation  distance,  approximately.  200-300  meters.  Table  1 
lists  our  estimates  of  the  limitations  to  angular  resolution  imposed  at 
VHF-UHF  frequencies  by  the  auroral  ionosphere. 

Table  1 

Limitations  to  Angular  Resolution  Imposed  by  Auroral  Ionosphere 

Frequency  in  Mcs  50  100  200  300 

Minimum  Angular  Separation  in  Radians  .12  .03  .0075  .003 

2.  3  Limitations  of  Angular  Resolution  Imposed  by  the  Troposphere 

At  the  time  of  the  initiation  of  the  Wide  Band  program,  there  was 
a  definite  paucity  of  available  data  on  the  magnitude  of  tropospheric  path 
fluctuations.  In  the  first  Wide  Band  report,  extensive  use  was  made  of 
the  path  length  fluctuation  data  obtained  by  the  National  Bureau  of  Standards 
from  «  10  mile  path  in  Colorado.  Since  that  time  NBS  has  collected 
additional  Colorado  data  and  extensive  Hawaiian  data.  Figure  1  is  a 
composite  of  the  NBS  data  reported  by  Thompson.  Janes  and  Kirkpatrick  (I960). 
The  data  presented  are  median  spectral  data  and  the  variance  of  the 
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Figure  I  SPECTRUM  OF  APPAREHT  PATH  LENGTH  FLUCTUATIONS  IN  THE  TROi'»OSPHERE 

(INOHPSON.  JANES,  KIRKPATRICK) 
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spectral  density  estimates  about  this  median  value  might  very  well 
be  approacimately  10  decibels.  The  important  points  to  note  in  these 
data  are  the  approximate  variation  of  the  spectral  density  with  the 
inverse  2. 6  power  of  the  fluctuation  frequency;  the  fact  that  the  median 
Hawaii  data  are  approximately  10  db  greater  than  Uie  median  Colorado 
data  and  finally  that  the  spectral  density  data  pertaining  to  the  passage 
of  a  frontal  system  at  the  Colorado  site  are  20  db  greater  than  the 
median  of  the  Colorado  data. 


In  order  to  indicate  reasonable  limits  to  the  spectral  density 
of  path  length  fluctuations  we  have  drawn  the  d:>Hed  line  in  Figure  1 
which  encompasses  all  of  the  median  experimental  data  and  should  be 
representative  of  the  majority  of  the  experimental  runs.  In  Figure  Z, 
we  present  the  expected  value  of  the  spectrum  based  upon  a  five  second 
data  sample  (see  Appendix  A).  The  area  under  this  curve,  the  variance 
of  the  path  length  fluctuations  is  11.4  x  10**^  (parts /million)^  If  we 
adopt  as  a  model  of  the  troposphere  a  statistically  uniform  medium 
extending  from  ground  up  to  ten  kilometers  altitude  and  a  lower 
limit  to  the  elevation  angle  of  six  degrees,  then  the  maximum  tropospheric 
path  length  will  be  less  than  100  kilometers  long.  Consequently  we 
expect  that  the  root  mean  square  phase  fluctuations  on  a  |ong  tropospheric 
path  will  be  less  than 
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Figure  2  EXPECTED  VALUE  OF  SPECTRUM  OF  SINGLE  PATH  PHASE  FLUCTUATIONS 
BASED  ON  A  5  SECOND  DATA  SAMPLE 


FLUCTUATION  OF  AFPARBIT 
PATH  LENOTM  IN  (PARTS 
PER  MILLI0N)VUNIT  SANDWIOTH 


The  greater  the  frequency,  the  greater  will  be  the  rme  phase 
fluctuations.  At  X^band,  the  highest  frequency  considered  in  this  report 

•  JTTmo^’ 

V 

^  O-i?/ 

Since  our  previous  discussion  of  ionospheric  phase  fluctuations 

% 

at  X-band  indicates  that  there  are  less  than  1  square  radian,  we  do  not 

expect  resolution  to  be  impaired.  It  should  be  further  noted  that  even 

if  the  phase  fluctuations  were  larger  than  one  square  radian,  the  amount 

of  beam  spreading  (due  to  the  continuous  portion  of  the  angular  spectrum) 

might  be  quite  small.  The  correlation  distance  of  the  phase  fluctuations 

(derived  in  Appendix  A)  is  of  the  or  ier  of  30  meters.  At  X>band.  the 

continuous  portion  of  the  spectrum  i  '  therefore  C  3  ^  -  1  milliradian. 

d.  3vio' 

2.4  The  Wide  Band  Satellite  anf.  Aircraft  Interferometer  Experiments 
It  would  appear  therefore  that  tropospheric  and  ionospheric  effects 
should  not  limit  angular  resolution.  In  order  to  prove  these  conclusions 
under  actual  field  conditions  an  experiment  utilizing  satellite  transmissions 
at  2250  Mcs  was  proposed.  Interferometers  having  baselines  of  100  and 
1000  meters  in  North-South  and  East-West  directions  were  installed  at 
a  field  site  near  Buffalo.  New  York.  The  equipment  is  fully  described 
in  Wide  Band  Report  No.  UB-1363-P-101,  The  satellite  transmitter  was 
designed  with  both  the  interferometer  experiment  and  wide  bandwidth 
transmission  e  .periment  (to  be  described  in  Section  4  on  Range  Resolution) 
in  mind. 
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The  interferometer  experiment  had  as  its  goal  the  measurement 


f 


of  the  phase  fluctuations  incurred  on  a  combined  tropospheric  and 
auroral-ionospheric  path.  In  order  to  accomplish  the  auroral 
measurements^  a  polar  orbiting  satellite  was  required.  Our  initial 
plan  was  to  include  the  Wide  Band  transmitting  apparatus  on  an 
experimental  Midas  satellite.  By  the  time  it  became  apparent  that 
a  satellite  experiment  was  not  immediately  forthcoming,  the  installation 
of  ground  equipment  was  complete  and  a  system  checkout  program 
utilizing  aircraft  transmissions  was  initiated.  Although  the  prime 
purpose  of  the  aircraft  program  was  to  check  out  the  interferometer 
system,  it  was  recognized  that  this  program  would  also  offer  an 
opportunity  to  measure  the  tropospheric  phase  fluctuations  encountered 
on  a  rotating  line  of  sight  path  in  contrast  with  the  NBS  experiment 
already  mentioned,  which  employed  a  path  whose  coordinate  end 
points  were  fixed.  Detailed  results  of  this  program  are  presented 
in  Appendix  C  but  the  general  results  can  be  summarized  as  follows. 

On  a  propagation  path  which  extended  from  ground  level  to  approxi¬ 
mately  40,000  feet  and  whose  length  was  60  nautical  miles,  the  measured 
mean  squared  phase  fluctuations  were  very  closely  approximated 
by  the  phase  fluctuations  associated  with  the  uncertainty  of  actual 
aircraft  position.  The  difference  between  our  estimates  of  phase 
fluctuation  and  the  phase  fluctuations  associated  with  unknown  aircraft 
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position  were  far  less  than  the  estimates  of  the  worst  possible 
phase  fluctuation  conditions  illustrated  in  Figure  2.  We  conclude 
therefore  that  Figure  2  probably  represents  a  pessimistic  estimate 
(which  it  was  designed  to  do)  of  the  phase  fluctuations  associated  with 
the  troposphere  and  therefore,  at  least  at  frequencies  equal  to  or 
less  than  10,000  Mcs,  the  troposphere  will  not  limit  angular  resolution. 


3.  VELOCITY  RESOLUTION 


The  ability  of  a  radar  system  to  resolve  two  targets  on  the 
basis  of  differing  Doppler  shifts  is  clearly  a  function  of  the  purity 
of  the  transmitted  waveform*  and  the  duration  of  the  observation.  It 
la  equally  a  function  of  the  stability  of  the  propagation  path.  Consider 
a  monochromatic  plane  wave  reflected  from  a  target  whose  radial 
velocity  is  .  In  the  absence  of  phase  fluctuations  due  to  the  propagation 
path,  the  returned  signal  is  of  the  form 

where 


If  there  are  random  phase  fluctuations,  the  received  signal  is  of  form 


tii)  = 


where  '^r;.  )  is  the  random  phase  fluctuation  associated  with  the 

propagation  path.  The  spectrum  of  the  received  signal  can  be  obtained 
from  the  Fourier  transform  of  the  autocorrelation  function  ^ 


,C 


f 
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where  defined  and  C0|(  i« 

the  Doppler  shifted  frequency.  As  before,  if  ^(()  normally 

distributed,  the  density  functim  of  ^(j)^  ('{)  !•  given  by  p 

Hr 


p(^4')  - 


f  -r  ■  j 


where  again  is  the  variance  of  ^^./Wand  »• 

normalized  autocorrelation  function  of  (t^) 


the 


/2r  ( T)  (^d.  0-7 frj ]j 


And  again,  in  a  manner  identical  to  that  indicated  in  Appendix  A,  the 
spectrum  of  is 


(r) 


cos  CujrcsscoT' 


Propagation  phase  path  fluctuations  ttusract  power  from  the 
Doppler  shifted  line  in  the  received  power  s^  ectrum  and  spread  it  into 
a  continuous  distribution  centered  around  the  Doppler  shifted  line.  So 
long  as  the  mean  squared  phase  fluctuations  are  small  (less  than  one  square 
radian)  the  power  in  the  "smeared"  continuous  spectrum  will  be  much 
less  than  the  p'wer  in  the  Doppler  shifted  line  and  velocity  resolution 
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of  two  target!  of  equal  or  nearly  equal  crosi  section  can  be  effected 
to  the  extent  permitted  by  signal-to*noise  ratio,  integration  time  and 
other  system  parameters. 

Consider  that  we  desire  to  resolve  targets  whose  velocities 
differ  by  0. 1  feet  per  second  (.03  meters /second).  The  differential 
Doppler  shift  will  be 

z  ( A./r)  _  ,  ot 

*  ~Tr  ~ 

Table  11  lists  the  differential  Doppler  shift  as  a  function  of  operating 
wavelength.  Also  listed  is  the  miniirum  time  necessary  for  any 
system  to  resolve  the  indicated  frequency  difference. 


Table  11 


Operating  Wavelength  in  Centimeters 

30 

10 

3 

Frequency  Difference  cps 

0.2 

0.6 

2.0 

Minimum  Observation  Time  seconds 

10 

3.3 

1 

For  an  averaging  time  of  5  seconds,  we  have  already  seen 
that  the  mean  squared  phase  fluctuations  atX-band  are  less  than  0.394 
square  radians.  Consequently  we  do  not  expect  that  tropospheric  phase  fluctu> 
ations  alone  will  limit  the  ability  to  achieve  velocity  resolution  of  0. 1 
feet  per  second. 
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4.  RANCfE  RESOLUTION 


Range  resolution  is  the  ability  of  a  radar  system  to  distinguish 
between  two  targets  of  equal  cross  section  which  differ  only  in  the  range 
coordinates.  If  the  two  targets  are  point  sources  and  have  identical 
radar  cross  sections,  then  the  range  resolution  achievable  by  a  radar 
system  is  slightly  in  excess  of  (1/2)  where  c  s  velocity  of  propagi^tion 
and  ^  is  the  duration  of  the  radar  pulse. 

Radar  range  resolution  is  a  function  of  the  bandwidth  of  the 
radar  system.  Increased  range  resolution  is  achieved  by  radiating 
a  narrower  pulse.  The  narrower  the  radiated  pulse,  the  broader  is 
the  radiated  spectrum  and  hence,  the  bandwidth  of  the  radar  system. 

« 

The  limits  to  the  achievable  range  resolution  are  the  limits 
imposed  upon  the  usable  bandwidth.  The  real  atmosphere  of  the  earth 
will  impose  bandwidth  limitations  because  the  atmosphere  is  both 
dispersive  and  non-uniform.  These  considerations  will  be  examined  in 
order. 

4. 1  The  Dispersive  Atmosphere 

The  presence  of  free  electrons  in  the  ionosphere  (a  region 
extending  from  approximately  60  kilometers  in  altitude  upwards  to  the 
limits  of  the  earth's  atmosphere)  modifies  the  refractive  index  such 
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60 


where 


M  s  refractive  index 

hj  -  lo'cal  electron  density  (number  per  cubic  meter) 

•3 

s  mass  of  an  electron  =  9  x  10**  kilograms 

_o 

^  s  charge  of  an  electron  s  1.  59  x  10  coulombs 
40  -  angular  radio  frequency 

-12 

^  =  capacitivity  of  free  space  =  8.854  x  ID  farads/meter 


It  is  to  be  noted  that  the  refracti\re  index  hi  the  ionosphere 
is  a  function  of  frequency.  Consequently  each  of  the  waves  (of  differing 
frequency)  which  comprise  the  radio  frequency  pulse  will  travel  with 
a  different  velocity.  The  pulse,  after  travelling  through  the  dispersive 
medium,  will  be  distorted.  The  distortion  will  be  evidenced  by 
an  increase  of  apparent  pulse  length  as  well  as  a  change  in  the  general 
shape  of  the  pulse. 


22 


In  the  Phase  1  Report  of  Project  Wide  Band  (CAL  Report 
No.  UB-1363-P-3,  March  I960)  we  reported  on  the  results  of  the 
effects  of  dispersion  upon  pulse  shape  as  a  Rinction  of  carrier  frequency, 
percentage  bandwidth  and  integrated  electron  density  along  the 
propagation  path.  Asa  result  of  this  investigation  we  concluded  that 
if  dispersion  alone  were  the  limiting  factor  to  pulse  shape,  then  (for 
frequencies  in  excess  of  4000  Mcs)  bandwidths  of  500  Mcs  could  be 
achieved  without  requiring  corrective  action.  It  also  appeared  that  if 
corrective  filters  were  employed,  500  Mcs  bandwidth  propagation 
could  be  accomplished  at  a  carrier  frequency  as  low  as  2000  Mcs. 

4.  2  Non-Uniform  Atmosphere 

The  fact  that  the  refractive  index  in  the  ionosphere  as  well 
as  in  the  troposphere  is  non-uniform  may  limit  the  range  resolution 
at  certain  frequencies.  To  the  extent  that  multipath  propagation 
(discussed  in  the  Section  on  Angular  Resolution)  is  important,  it 
may  limit  the  achievable  range  resolution  since  energy  travelling 
over  non  line*  of- sight  paths  will  arrive  at  times  later  than  the  energy 
travelling  over  the  direct  path.  In  the  Phase  1  Report  on 
Project  Wide  Band,  we  indicated  that  we  were  not  certain  of  the 
multipath  situation  in  the  auroral  ionosphere.  As  a  result  of  the  Radio 
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Astronomy  Experiment  reported  in  the  first  Radio  Astronomy  Report 
CAL  Report  No.  UB-1363-P-4]  and  with  some  additional  detail  in 
Appendix  B,we  have  concluded  that  range  resolution  o*  i  foot  (a  bandwidth 
of  500  Mcs)  is  possible  in  the  presence  of  the  auroral  ionosphere  if 
the  carrier  frequency  is  well  in  excess  of  500  Mcs. 

It  should  be  emphasized  that  the  above  conclusion  is  based 
upon  a  defmition  of  resolution  which  implies  that  the  target  radar 
cross  sections  are  identical  and  that  the  phase  fluctuation  statistics 
are  not  strongly  frequency  dependent.  We  shall  discuss  the  first 
limitation  in  the  next  section  of  this  report  and  shall  indicate  the  extent 
of  the  resolution  problem  for  unequal  targets  in  greater  detail. 

During  the  course  of  this  discussion,  the  effects  of  the  frequency  dependence 
of  the  phase  fluctuation  statistics  will  become  more  clear. 


5.  RESOLUTION  OF  TARGETS  OF 
UNEQUAL  RADAR  CROSS  SECTION 

In  the  previous  discussions  we  have  been  able  to  give  our 
estimates  of  the  resolution  lin.its  imposed  by  the  atmosphere  in 
angle,  range  and  velocity  space  when  resolution  is  defined  in  terms 
of  targets  of  equal  cross  sections.  In  recent  years,  the  specifications 
for  more  sophisticated  radar  systems  l^ave  pre* supposed  the  existence 
of  systems  which  can  resolve  targets  whose  radar  cross  section  may 
differ  by  more  than  30  db.  Thfj  effects  of  the  atmosphere  upon  the 
resolution  capabilities  of  such  systems  poses  a  problem  of  quite 
different  magnitude  than  for  the  types  of  system  which  have  previously 
been  considered.  A  simple  example  might  serve  to  illustrate  the  poiht. 

Consider  the  angular  resolution  problem  treated  in  Section  2. 

We  have  seen  that  when  the  mean  squared  phase  fluctuations  are  CT  , 
power  is  abstracted  from  the  wave  arising  from  the  target  direction 
and  is  spread  over  a  range  of  angles  of  roughly  It  • 

The  available  data  indicate  that  the  expected  value  of  the  mean  squared 
phase  fluctuations  is  less  than  0.  394  radians  and  consequently  the  full 
width  of  the  continuous  spectrum  is  approximately  one  milliradian  at 
X-band.  The  power  associated  with  the  continuous  spectrum  is  (  /**  ^  ) 

so  that  for  example  the  power  in  a  range  of  angles  ^0  ,  of  approximately 

1/10  milliradians  is  .032  units.  The  situation  is  depicted  in  Figure  3 
where  the  resultant  undiffracted  wave  is  13  db  more  intense  than  the 
power  in  a  tenth  milliradian  beamwidth.  Consequently  it  would 


Figure  3  ANGULAR  SPECTRUM  AT  X>BAND 

ANTENNA  BEAHWiDTH  »  0. I  MILL  I  RADIANS 
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be  very  difficult  to  resolve  two  targets,  separated  by  a  one-half 
milliradian  whose  radar  cross  sections  differed  by  more  than  13  db. 

The  above  results  are  dependent  upon  the  mean  squared  phase  fluctuations 
being  approximately  0.394  square  radians,  the  correlation  distance 
being  30  meters,  and  the  antenna  beamwidth  being  0. 1  milliradians. 

Another  example  may  illustrate  the  dependence  of  the  resolution 
criterion  upon  these  parameters.  If  all  values  of  the  parameters  are 
the  same  as  before  except  that  the  correlation  distance  of  the  atmospheric 
phase  fluctuations  is  increased  to  100  meters,  then  the  power  density 
per  unit  0. 1  milliradian  beamwidth  will  be  three  times  the  previous 
value  and  the  width  of  the  continuous  spectrum  will  be  one -third  of  a 
milliradian.  Under  these  circumstances,  a  one-tenth  milliradian  system 
will  have  difficulty  in  resolving  two  targets,  separated  by  0.16  milliradians, 
when  their  cross  sections  differ  by  more  than  15  db.  If  the  averaging 
time  required  to  make  the  measurements  is  less  than  the  5  seconds 
we  have  used  in  this  report,  the  expected  variance  of  the  phase  fluctuations 
will  be  reduced  and  the  resolution  capabilities  of  the  system  will  be 
improved.  It  should  also  be  noted  that  if  antenna  size  presents  no 
insuperable  engineering  problems,  decreasing  the  operating  frequency 
while  increasing  the  antenna  aperture  will  reduce  (by  the  frequency  ratio 
squared)  the  variance  of  the  phase  fluctuations  and  hence  improve  the 
angular  resolution  of  a  system. 
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Similar  considerations  hold  for  the  velocity  resolution 
capabilities  of  a  radar  system.  In  this  case  the  "correlation  time" 
of  the  phase  fluctuations  takes  the  place  of  the  "correlation  distance" 
used  in  the  previous  angular  resolution  discussion.  Reducing  the 
averaging  time  or  the  frequency  of  operation  may  improve  the  resolution 
capabilities  of  proposed  systems. 

The  ability  of  a  system  to  resolve  two  targets  of  vastly  unequal 
cross  section  which  differ  only  in  range  is  a  function  of  what  has  come 
to  be  called  "range  sidelobes".  Range  resolution  of  approximately 
one  foot  (500  Mcs  bandwidths)  is  usually  only  considered  for  pulse 
compression  systems.  In  order  to  achieve  resolution  of  this  order  a 
swept  frequency  pulse  approximately  5  or  10  microseconds  is  collapsed 
down  to  the  desired  pulse  length.  Rather  precise  shaping  of  the  amplitude 
and  phase  characteristics  of  this  pulse  is  required  to  achieve  range 
sidelobe  levels  of  perhaps  25-30  db.  Strong  signal-to-noise  ratios  are 
required  if  the  phase  of  the  signal  is  to  be  a  meaningful  quantity.  If 
the  phase  path  fluctuations  are  exactly  governed  by  a  relationship  of 
the  sort 
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where  ^ U  the  inetenteneoae  pheee  ehift  et  time  ^  end  wevelepgthX  | 

A^tl)  it  the  random  phase  path  fluctuation  and  ia  independent 
of  frequency 

then  the  only  way  in  which  phase  path  fluctuations  can  influence  the 

shape  of  the  resultant  collapsed  pulse  ia  to  have  the  random  phase 

/> 

path  fluctuation  change  in  the  time  period  u  of  the  uncollapsed  pulse. 

It  should  be  further  noted  that  a  linear  change  in  phase  path  with  time 
may  cause  a  range  error  but  cannot  cause  an  increase  in  range  sidelobes. 

The  available  propagation  data  indicate  that  the  changes  of 
phase  path  fluctuation  in  time  period  less  than  10  microseconds 
are  very  ^mall  and  indeed  are  beyond  the  measurement  capabilities 
of  available  equipment.  We  are  led  to  believe  then  that  if  the  uncollapsed 
pulse  is  less  than  10  microseconds  long,  with  the  signal«to-noise  ratios  and 
linear  dynamic  range  of  available  wide  band  radar  equipment,  changes  in 
phase  path  fluctuations  during  the  pulse  time  is  of  minor  importance. 

What  may  be  of  greater  importance  is  the  fact  that  the  path 
length  fluctuations  at  two  widely  separated  frequencies  may  be  uncorrelated 
and  cannot  therefore  be  represented  by  (5.1).  Two  possible  mechanisms 
come  readily  to  mind  which  may  lead  to  this  situation.  The  first  mechanism 
involves  the  dispersive  properties  cf  the  ionosphere  which  could  cause 
two  frequencies,  separated  by  500  Mcqto  travel  two  slightly  different 
ray  paths  and  therefore  encounter  two  differing  time  histories  of  path 
fluctuation. 
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One  might  expect  that  the  contribution  of  the  very  largest  irregularities 
of  refractive  index  might  be  common  to  both  radio  frequencies  but 
that  path  length  fluctuations  of  the  smaller  irregiilarities  would  be 
uncorrelated. 

The  second  mechanism  has  to  do  with  the  scattering  properties 
of  the  irregularities  of  refractive  index.  If  there  is  a  cutoff  size  in 
the  spectrum  of  irregularities,  the  smallest  irregularities  will  scatter 
the  lower  of  the  two  frequencies  over  a  slightly  wider  angle  than  they 
will  the  higher  of  the  two  frequencies.  Since  the  resultant  phase  at  a 
point  on  the  ground  is  the  sum  of  all  such  forward  scattered  wavelets, 
the  two  path  length  histories  may  be  slightly  different  at  the  two  different 
frequencies. 

Both  of  the  suggested  mechanisms  for  uncorrelated  path  length 
fluctuations  recognize  that  such  uncorrelated  path  length  fluctuations 
will  be  primarily  due  to  the  smaller  refractive  index  irregularities.  On»» 
would  therefore  expect  that  the  higher  frequency  components  of  the 
phase  fluctuations  would  be  uncorrelated.  The  National  Bureau  of 
Standards  has  completed  a  series  of  measurements  on  a  47  kilometer 
line-of- sight  over-water  path  at  Eleuthera  Island  in  the  Bahamas  using 
frequencies  of  9400  and  9200  Mcs.  (Janes  et  al.  1963).  Although  the 
presence  of  the  sea  surface  as  a  reflector  tends  to  obscure  the  applicability 
of  the  results,  the  data  show  that  the  phase  history  of  the  two  signals 
is  uncorrelated  at  spectral  frequencies  greater  than  one  cycle  per  second. 
This  is  the  only  experiment  instrumented  to  measure  uncorrelated  path 
length  fluctuations  as  a  function  of  frequency  separation.  Such  path  length 
fluctuations  may  set  a  lower  limit  of  range  sidelobe  suppression. 
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5, 1  The  Wide  Band  Satellite  Experiment 

In  order  to  measure  the  effects  of  uncorrelated  phase  path 
fluctuations  for  frequencies  differing  by  500  Mcs  on  auroral  ionospheric 
and  long  tropospheric  paths  a  satellite  experiment  was  proposed  and 
initiated.  The  basic  experiment  consisted  of  sending  a  carrier  and  3  pairs 
of  sidebands  separated  by  100«  200,  and  300  Mcs.  These  signals  were 
to  be  received  on  the  ground  and  the  sidebands  were  to  be  separately 
mixed  with  the  carrier  to  obtain  pairs  of  frequencies  of  50,  100,  and 
150  Mcs.  The  phase  difference  between  these  pairs  of  frequencies 
was  to  be  recorded  and  analyzed.  The  dispersive 

properties  of  the  atmosphere  would  cause  the  recorded  phase  difference 
to  vary  slowly  and  regularly  (the  phase  difference  is  proportional  to  the 
integrated  electron  density  along  the  propagation  path).  Measurements 
of  the  phase  fluctuations  about  the  regular  curve  should  lead  to  estimates 
of  the  mean  squared  phase  fluctuations  as  a  function  of  frequency  separation. 
In  this  experiment  the  difference  phase  fluctuation  variance  would  be 
measured  at  only  three  frequency  separations. 

We  believed  that  it  was  most  important  to  make  these  measure¬ 
ments  through  the  auroral  ionosphere  so  that  the  satellite  vehicle  should  be 
in  a  near-polar  orbit.  The  carrier  frequency  of  the  Wide  Band  Satellite 
transmitter  was  also  to  be  used  for  the  satellite  interferometer  experiment 
described  in  Section  2  of  this  report. 


31 


The  choice  of  carrier  frequency  and  actual  radiated  bandwidth 
of  the  satellite -borne  transmitter  was  determined  by  the  availability 
of  hardware  which  could  operate  in  the  satellite  environment.  A  study 
and  evaluation  of  the  available  designs  of  wideband  transmitters  was 
presented  in  the  Phase  1  Wide  Band  report.  On  the  basis  of  the 
information  then  available  a  frequency  modulated,  voltage  tuned 
magnetron  was  selected  for  the  transmitter.  A  subcontract  was  let 
with  Radiation,  Inc.  of  Melbourne,  Florida  to  furnish  this  transmitter  in 
conformance  with  specifications  detailed  in  the  Phase  1  Wide  Bana 
Report.  During  the  course  of  this  development  program  difficulties 
with  magnetron  performance  were  experienced  and  the  subcontract 
transmitter  program  was  discontinued  in  favor  of  a  new  design  making 
use  of  TWT's  that  had  been  developed  for  satellite  use  since  initiation  of 
this  project. 

When  at  a  later  time  it  became  apparent  that  the  Wide  Band  Satellite 
Transmitter  could  not  be  acconrmodated  "piggyback"  on  existing  Air  Force 
satellite  programs,  attention  was  directed  towards  the  design  of  a 
complete  satellite  package  including  solar  cells,  command  and  control, 
telemetry,  and  atitude  stabilization.  The  results  of  this  design  study 
are  reported  in  CAL  Report  UB-1363-P-Z01  entitled  "Spacecraft  for  Radar 
Propagation  Experiments-Design  Study  Report"  dated  30  June  1962. 

The  spacecraft  design  was  adaptable  to  larger  launching  vehicles  but 
it  appeared  that  the  Air  Force  Blue  Scout  launch  vehicle  would  suffice. 
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The  funds  required  for  the  spacecraft  development  program  were 
unavailable  and  as  a  result,  no  active  satellite  experiments,  have  been 
performed. 

5,  2  Radar  Studies  of  Passive  Satellites 

The  original  Wide  Band  contract  work  statement  called  for  an 
investigation  of  the  feasibility  of  establishing  the  atmospheric  limits  on 
resolution  from  radar  measurements  of  orbiting  satellite  vehicles. 

At  that  time  (I960)  it  was  apparent  that  existing  radar  systems  having 
sufficient  sensitivity  to  measure  satellite  echo  strength  as  a  function 
of  time  did  not  have  the  necessary  bandwidth  to  enable  one  to  estimate 
range  resolution  limitations.  The  cost  of  installing  coherent  transmitting 
and  receiving  systems  which  could  utilize  satellite  reflections  at  say  three 
frequencies  separated  by  approximately  300  to  500  Mcs  was  deemed 
inordinately  high.  Consequently  attention  was  focussed  on  one  way 
propagation  measurements  using  a  satellite -borne  transmitter. 


5,  3  Additional  Areas  of  Investigation 


During  the  course  of  the  Wide  Band  program  we  were  directed 
to  analyze  the  atmospheric  propagation  effects  upon  a  dispersed  antenna 
array  •  satellite  communication  system.  The  final  report  for  this 
preliminary  investigation  is  submitted  as  Appendix  D  of  this  report. 
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6.  SUMMARY 


The  limitations  imposed  upon  radar  resolution  as  a  consequence 
of  propagation  through  the  atmosphere  have  been  defined  in  terms 
of  the  ability  of  a  radar  system  to  discriminate  between  two  targets  of 
equal  radar  cross  section.  In  so  doing  we  have  attained  the  basic 
objectives  of  Project  Wide  Band.  The  experimental  data  obtained  under 
the  radio  astronomy  portion  of  the  Wide  Band  program  have  been  most 
valuable  in  defining  the  limits  to  which  the  auroral  ionosphere  will  limit 
resolution  as  a  function  of  frequency.  The  recent  data  on  tropospheric 
fluctuations, obtained  by  the  National  Bureau  of  Standards  under  a  variety 
of  weather  and  climatic  conditions,  provides  the  basis  for  estimates  of  the 
limits  of  resolution  imposed  by  tropospheric  propagation.  The  fact 
that  our  estimates  are  pessimistic  and  are  designed  to  reflect  the 
worst  propagation  conditions  has  been  confirmed  in  part  by  our  wide 
baseline  interferometer -aircraft  experiment. 

The  limits  to  radar  resolution  in  the  case  of  two  targets  of 
unequal  cross  section  cannot  be  simply  defined  since  the  end  results 
are  dependent  upon  the  individual  radar  system  time  constants  and 
parameters.  However,  the  method  of  evaluating  the  limitations  of 
augular  and  velocity  resolution  (for  unequal  targets)  of  a  system  has 
been  demonstrated  through  an  example  in  Section  5.  We  also  showed 
that  to  the  ext^int  that  propagation  path  length  fluctuations  are  independent 
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of  frequency,  range  resolution  will  not  be  limited  by  atmospheric  effects. 
However,  we  noted  that  at  least  part  of  the  path  length  fluctuations 
at  two  separated  frequencies  will  be  uncorrelated.  If  the  variance  of 
these  uncorrelated  phase  fluctuations  is  small  enough,  the  attendant 
increase  in  "range  sidelobe  levels"  may  not  be  noticeable  in  the  sidelobe 
levels  which  will  exist  in  an  unperturbed  system.  Since  there  exist  no 
measurements  of  the  uncorrelated  phase  fluctuations  along  a  complete 
propagation  (auroral- ionospheric  plus  tropospheric)  path  we  have  no  infor¬ 
mation  on  which  to  base  an  estimate  of  the  resolution  limitations  which 
might  be  imposed.  We  feel,  however,  that  for  a  1-foot  pulse -length 
system  operating  at  frequencies  of  4000  Mcs  or  greater,  the  resolution 
limitations  will  only  be  associated  with  targets  which  have  grossly 
disparate  radar  cross  sections,  say  25-30  db. 

We  conclude  with  a  list  of  resolution  capabilities  for  targets 
of  equal  cross  section  which  we  feel  can  be  obtained  under  the  worst 
propagation  conditions  associated  with  the  earth's  atmosphere.  To  the 
extent  that  the  auroral  ionosphere  is  representative  of  the  conditions 
to  be  expected  as  a  result  of  a  high  altitude  nuclear  detonation,  then 
these  estimates  will  hold  for  that  condition, 

1.  Angular  Resolution  -  Angular  resolution  of  0,1  milliradians, 
at  frequencies  greater  than  1000  Mcs  can  be  attained  for  targets  of  equal 
cross  section.  Our  estimates  are  that  this  degree  of  resolution  can  be 
obtained  at  X-band  when  the  target  cross  sections  differ  by  as  much  as  13  db. 
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Z,  Velocity  Resolution  -  Velocity  resolution  of  0. 1  feet  per 
second  can  be  attained  for  targets  of  equal  cross  section.  We  estimate 
that  this  degree  of  resolution  can  be  obtained  at  X-band  when  the 
target  cross  sections  differ  by  as  much  as  13  db. 

3.  Range  Resolution  •  Rahge  resolution  of  1.0  foot  at  frequencies 
greater  than  4000  Mcs  can  be  attained  for  targets  of  equal  cross  section. 
The  limitations  to  resolution  for  targets  differing  greatly  in  radar 
cross  section  cannot  be  well  defined  but  may  be  equal  to  1.0  foot 
for  targets  whose  cross  sections  differ  by  as  much  as  30  db. 

The  above  estimates  are  based  upon  essentially  infinite  signal- 
to-noise  ratio.  The  range  and  velocity  resolution  achievable  with 
practical  radar  systems  may  very  well  be  limited  by  the  limits  to  signal* 
to -noise  ratio  and  the  linear  dynamic  range  of  operation  inherent  in 
any  practical  radar  system. 
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Appendix  A 

Mathematical  Derivation* 


Pha*e  Changing  Diffraction  Screen 

Consider  that  the  electric  field  inteneity,  meaeure  along  a 
plane,  of  a  radio  wave  ie  given  by 

'  <ros I /vf  ^ 

and  that  the  field  measured  at  a  point  X  i  AK  is 

o 


=  cos[ cjrf  (()] 


The  spatial  autocorrelation  function  of  the  electric  field  as  a 


function  of  y  ,  the  distance  between  the  measuring  points 


C'f) 


If  the  are  normally  distributed,  it  can  be  shown  that 

the  density  function  of  the  difference  e  ^  ^  is  normally 


distributed  with 


where 


O' 


is  the  probability  density  function  of 

is  the  variance  of  the  phase  fluctuations 

is  the  normalized  spatial  autocorrelation  function  of 
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The  angular  power  spectrum  is  the  Fourier  transform  of  the  spatial 
correlation  function  of  the  electric  field 

/ Ffs)l S(s)  i  A'  'CS) ] 

where  >  ^  is  the  sine  of  the  angle  of  arrival 

In  the  case  considered  by  Bramley 

AW*  ^/’/-f?jy 

where  /q  is  the  "correlation  distance"  measured  in  wavelengths. 


so  that 
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When  (J  /  »  terms  in  the  series  for  ^  out  to  approximately  t)*- ^  (T 
must  be  taken.  The  half  width  of  the  continuous  iipectrum  f-^CS )  will 


then  be 
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T  yo  ad 
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a  random  variable 


Given  f(^)^  ^  Cosj  CDt»'6 
Vis  j 

we  will  define  the  autocorrelation  function  of  a  finite  eample  length  7^  as 
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If  we  aek  for  the  exacted  value  of  the  correlation  function  and  recogniaing 


that  in  the  equations  the  values  of 


and  j  are  equally  probaUe 


escpeeted.  value  of 


a 


The  variance  is  seen  to  be  composed  of  the  ^ower  in  each  frequency  band 
which  would  exist  if  we  had  an  infinite  record  length  ^  ,  but 

with  each  power  density  weighted  by  a  filter  function 


jfa).-  /-  r 
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For  high  frequencies* 


is  very  smell  end  the 


J 


high  frequency  spectre!  densities  contribute  their  full  weight  to  the 
totel  verience.  Computing  the  veriance  from  e  limited  time  average  is 
equivalent  to  passing  the  data  through  a  high  pass  filter  whose  3  db 
cutoff  occurs  when 

^  "  fe.  ^ .  O-  7o7 
rrfj 


yfc  -  ^  - 


For  frequencies  below  and  up  to  the  cutoff  frequency 
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so  that  the  filter  function  behaves  as 
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Appendix  B 

Radio  Star  Fadeouts  and  Radio  Propagation 


Interferonneter  records  of  intense  radio  sources  have  sometimes 
shown  the  loss  of  the  interference  patterns  for  periods  of  time  ranging 
from  a  few  minutes  to  several  hours.  Tnis  phenomenon,  now  called  a 
radio  star  "fadeout,  "  was  first  reported  by  Huntley  (1953)  and  later  by 
Koster  (1956,  1958  and  1963)  and  was  thought  to  be  characteristic 
primarily  of  the  equatorial  ionosphere.  More  recently,  similar  radio 
star  "fadeouts”  have  been  reported  by  investigators  viewing  radio  sources 
through  the  ionosphere  in  high  northern  magnetic  latitudes  (Little  (1959). 
Benson  (I960),  Flood  (1962),  Moorcroft  (1963),  Flood  (1963)).  Radio 
star  "fadeouts"  have  also  been  observed  in  Hawaii  (to  be  discussed  in  this 
paper),  and  severe  amplitude  scintillations  were  also  observed  by  a  400  Mcs 
pencil  beam  radiometer  located  in  Trinidad,  B.W.l. 

In  this  Giscussion  we  review  the  results  reported  in  Flood  (1963), 
and  present  new  results  obtained  from: 

1.  A  further  analysis  of  the  data  collected  during  i960  -  1961 

2.  An  analysis  of  radio  star  fadeouts  observed  in  Hawaii 

3.  A  period  of  intense  scintillation  of  the  radio  sources 
in  Cygnus  and  Cassiopeia  as  measured  by  a  400  Mcs 
radiometer  in  Trinidad,  B.W.l. 


Review  of  Previous  Results 


Phase  switched  interferometers  operating  at  52  and  147  Mcs  using 
baselines  of  300,  600  and  900  feet  East-West  and  1000  and  2000  feet  North 
South  were  employed  to  observe  the  radio  source  in  Cassiopeia.  At  the 
latitude  of  Buffalo,  New  York  (43”),  the  radio  source  is  circumpolar  and 
during  the  period  September  I960  to  September  1961,  24-hour  recordings 
of  the  interference  pattern  were  obtained  at  62  Mcs.  The  147  Mcs  equip¬ 
ment  was  initially  directed  so  that  Cass  A  was  monitored  around  the 
hours  of  lower  culmination.  After  April  1961,  the  antennas  were  modified 
so  that  24-hour  recordings  were  obtained  at  147  Mcs  as  well. 


1 
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In  Flood  (1963),  the  results  of  the  East- West  baselines  were  analyzed. 
Radio  star  fadeouts  were  shown  to  be  fairly  common  occurrences.  For 
reference  purposes.  Figure B>lis  an  example  of  the  fringe  pattern  at  52  Mcs, 
taken  during  a  "quiet"  uppei  culmination  period.  An  example  of  a  radio 
star  fadeout  is  given  in  Figure  B-2.  The  fadeout  is  indicated  by  an  arrow  in 
the  figure.  Notice  that  the  reduction  in  amplitude  of  the  interference  pattern 
is  a  function  of  the  antenna  baseline  -  the  reduction  in  amplitude  increases 
with  increasing  baseline  distance.  Absorption  therefore  cannot  be  a  prime 
factor  in  producing  a  fade.  Figure  ^3  isa  plot  of  the  number  of  occurrences 
of  fadeouts  as  a  function  of  .local  time.  Star  fadeouts,  while  occurring 
uniformly  thru  all  hours  of  local  time,  do  show  a  distinct  tendency  to 
maximize  in  the  evening  hours. 

Figure B4  isa  plot  of  the  occurrence  of  radio  star  fadeouts  as  a  function 
of  local  time  for  twelve  successive  months.  The  average  time  of  lower 
culmination  for  each  month  is  indicated.  It  is  clear  that  while  star  fadeouts 
are  most  common  around  the  time  of  lower  culmination  (which  advances  two 
hours  in  local  time  for  each  month),  there  is  a  secondary  process  clearly 
visible  during  the  autumnal  equinoctial  periods.  To  better  illustrate  this 
point,  Figure B>5is  a  plot  of  the  number  of  radio  star  fadeouts  as  a  function 
of  local  time  during  the  months  oi  April  and  October.  During  April,  when 
the  time  of  lower  culmination  occurs  during  the  early  evening  hours  (2200 
EST),  the  curve  of  occurrence  of  star  fades  is  single-humped  and  centered 
around  the  time  of  lower  culmination.  7'he  October  data  on  the  other  hand 
show  two  clearly  defined  periods  of  fadeouts,  centered  around  the  time  of 
lower  culmination  (1000  EST)  and  the  early  evening  hour  (2300  hours). 

Figures  B-5&  B6  are  interpreted  as  indicating  the  existence  of  two 
independent  processes  producing  the  radio  star  fade.  The  first  is  a  zenith 
angle  dependence  such  that  star  fades  are  far  more  common  during  the 
hours  of  lower  culmination  when  the  path  lengths  through  a  scintillation 
causing  ionosphere  are  longest.  The  second  process  appears  to  be  a 
phenomenon  centered  in  local  time  in  the  early  evening  hours.  Since  the 
time  of  lower  culmination  progresses  steadily  through  all  hours  of  local 
time  in  the  course  of  a  year,  the  two  processes  converge  to  give  a  single 
temporal  peak  in  the  early  evening  hours  in  April  and  diverge  to  produce 
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Figure  B- I  QUIET  DAY  FRINGE  PAHERN 


FIGURE  B-2  INTERFEROMETER  FRINGE  PATTERN  INDICATING  A  RADIO  STAR  FADE 
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Figure  B-3  OCCURRENCE  OF  RADIO  STAR  FADES  OBSERVED  ON  900  FOOT  EAST-WEST 
BASELINE  AT  52  MCS  SEPTEMBER  20,  I960  TO  SEPTEMBER  20,  1962 
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Figure  B>5  DIURNAL  OCCURRENCE  RADIO  STAR  FADES  IN  OCTOBER  AND  APRIL 
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Figure  B>6  NUMBER  OF  HOURS  HAVING  RADIO  STAR  FADEOUTS  AS  A  FUNCTION  OF  STAR 
POSITION,  DATA  TAKEN  SEPTEMBER  20,  I960  TO  SEPTEMBER  30,  1961 
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the  biroodel  curve  of  October.  It  is  also  clear  that  the  majority  of  the 
fadeouts  are  produced  by  the  first  mechanism  **  associated  with  the  longest 
paths  thru  the  scintillation  producing  region.  This  fact  is  well  illustrated 
by  Figure  B^vhich  is  a  plot  of  the  number  of  star  fadeouts  as  a  function  of 
local  sidereal  time.  The  data  are  quite  well  fitted  hy  a  secant  dependence 
of  the  zenith  angle. 

Correlation  with  Other  Geophysical  Data 

Attempts  have  been  made  to  correlate  the  occurrence  of  radio  star 
fades  with  other  geophysical  parameters.  Before  discussing  these  results 
it  is  well  to  keep  in  mind  the  fact  that  the  F- region  penetration  point  for  the 
path  Buffalo  to  Cassiopeia  at  lower  culmination  is  at  latitude  52.  5*N  (1000  km 
north  of  Buffalo)  where  the  magnetic  dip  is  80*.  The  £-region  penetration 
point  at  lower  culmination  is  at  latitude  47*N  (450  km  north  of  Buffalo) 
where  the  magnetic  dip  angle  is  approximately  76*.  Magnetic  K  indices 
recorded  at  Agincourt,  Ontario  (approximately  100  km  north  of  Buffalo) 
provide  the  records  of  magnetic  activity  closest  to  the  £  and  F  region 
penetration  points.  The  ionosonde  closest  to  the  ionospheric  penetration 
point  is  located  in  Ottawa,  Ontario  (45.  4*  N,  75.  9* W),  appr'/ximately  400 
kilometers  northeast  of  Buffalo,  New  York.  Auroral  backscatter  radars, 
at  49.  7,  143.  5  and  226  Mcs  are  operated  from  Buffalo  under  a  program  of 
auroral  research.  The  most  common  range  of  auroral  echoes  at  Buffalo 
lies  between  500  and  800  km  north.  In  general,  therefore,  there  are  no 
geophysical  data  taken  precisely  at  the  ionospheric  penetration  points  for 
the  star  at  lower  culmination,  when  radio  star  fadeouts  are  most  probable. 

In  view  of  the  great  distance  between  Agincourt  and  the  £  and  F 
region  ionospheric  penetration  points,  we  attempted  to  correlate  the 
occurrence  of  radio  star  fadeouts  for  those  months  when  upper  culmination 
(for  which  the  penetration  points  are  closest  to  Buffalo)  occurred  during  the 
period  1800  -  2400  hours.  To  do  this  the  data  from  the  months  of  September, 
October,  November  and  December  I960  were  used.  During  these  months, 
the  times  of  lower  culmination  differ  sufficiently  from  the  evening  maxima 
times  so  as  to  be  resolvable  ->  the  monthly  diurnal  variation  curves  are  more 
or  less  "two-humped." 
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It  it  tempting  to  tpeculate  that  this  secondary  evening  maximum  it 
associated  with  auroral  activity  along  the  propagation  path.  If  we  can 
establish  a  relationship  between  K  index  at  Agincourt  and  the  presence  of 
backscattered  auroral  echoes  observed  at  Buffalo,  we  may  be  able  to  use 
the  magnetic  data  as  an  index  of  auroral  activity.  To  this  end.  film  records 
of  auroral  backscattered  echoes  at  49.  7  Met  were  examined  for  their 
correlation  with  Agincourt  K  figure.  The  results  are  plotted  in  Figure  B>7 
where,  as  can  be  seen,  the  probability  of  observing  an  aurora  at  50  Mes 
iner  sases  rapidly  for  K  values  greater  than  3.  There  is  only  one  data 
point  each  for  K  =  6  and  K  =  7  but  the  trend  of  the  curve  is  clear.  It  is  also 
interesting  to  note  that  backscattered  auroral  echoes  have  been  observed, 
although  infrequently,  when  the  Agincourt  K  index  was  0  and  1.  For  the 
purposes  of  this  paper  we  shall  assume  that  an  Agincourt  K  Index  of  four 
or  greater  implies  auroral  conditions. 

The  data  for  the  months  September,  October,  November  and  December 
1960  were  searched  for  fadeout  periods  during  the  hours  1800  to  2400.  The 
results  are  tabulated  in  Table  B*l. 


Table  B-I 

Fadeouts  during  1800*  i400  Local  Time 
Average  Agincourt  K 

Index  Non-fadeout  Average  K  Index  Total  Number  Number  of  Events 


Month 

Periods  Fadeout  Period 

of  Events 

with  K2.4 

Sept. 

1.97  3.1 

18 

6 

Oct. 

3.  0  3.  0 

15 

6 

Nov. 

3.0  6.5 

7 

6 

Dec. 

2.5  3.0 

4 

2 

Average  K 
Average  K 
Twenty  out 

for  entire  period  (1800-  2400  local  time) 
for  fadeout  periods 

of  44  events  have  K  '^4. 

=  2.6 

3.6 

The  average  value  of  the  K  index  ie  noticeably  higher  during  the  fadeout 
periods  than  it  is  during  non-fadeout  times.  Approximately  one  half  the  fadeout 
periods  may  be  associated  with  a  strong  probability  of  aurora  (K^4). 
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OB  LOSS  2  METERS 


Figure  ^7  PROIABILITY  OF  OBSERVING  SO  NCS  AURAL  ECHOES  FROM 
BUFFALO,  N.Y.  AS  A  FUNCTION  OF  A6INC0URT  X  INDEX 


Figure  B>8  DB  LOSS  AT  6  METERS  PLOTTED  AGAINST  DB  LOSS  AT  2  METERS 


A  substantial  number  of  auroral  backscatter  film  records  for  the  month 
of  August  1961  were  available  for  comparison  with  star  fadeout  periods  between 
1800  -  2400  hours.  Table  B1  lists  a  comparison  of  the  auroral  and  fadeout  data. 

Table  B-U 

A  Comparison  of  Aurora  Backscatter  and  Star  Fadeouts 

Max  K 


Date 

Time 

Fadeout 

Figure 

Auroral  Echoes 

Aug. 

1 

1844 

Yes 

3 

Yes 

Aug. 

2 

1824 

Yes 

3 

Yes 

Aug. 

3 

1800-1900 

Yes 

4 

Yes 

Aug. 

3 

1900-2000 

Yes 

5 

Yes 

Aug. 

4 

1800-2040 

No 

1 

No 

Aug. 

5 

1814 

Yes 

3 

No  . 

Aug. 

7 

1800-2246 

No 

2 

Yes  - 

Aug. 

8 

1800-2300 

No 

3 

No 

Aug. 

9 

1800-2117 

No 

2 

No 

Aug. 

10 

1800-2030 

No 

4 

Yes  (2000-2015) 

Aug. 

11 

1800-2323 

Yes  (1830- 

1930)  4 

Yes  (1814-1825) 

Aug. 

21 

1800-2214 

No 

1 

No 

Aug. 

22 

1800-2215 

No 

1 

No 

Aug. 

24 

1800-2318 

No 

3 

No 

Aug. 

25 

1800-2400 

No 

4 

No 

Aug, 

29 

1830 

Yes 

5 

Yes 

Aug. 

30 

1800-1900 

Yes 

4 

Yes 

Aug. 

30 

1900-2000 

Yes 

5 

Yes 

Aug. 

30 

2000-2100 

Yes 

5 

Yes  (Weak) 

Aug. 

30 

2100-2200 

No 

5 

Yes 

Aug. 

30 

2300-2100 

Yes 

4 

Yes 

We 

conclude  Ihetefore  that  there  is 

a  strong  corre 

lation  between  the 

presence 

of  aurora  ci.d  the 

occurrence  of 

a  star  fadeout. 

We  further  hypothe- 

size  that  the  -rly  evening  maximum  of  star  fades,  present  in  the  diurnal 
variation  cvirves  of  the  nionths  around  the  autumnal  equinox, is  correlated  with 
the  presence  of  aurora,  ‘  aese  con(  lions  are  in  substantial  agreement  with 
those  of  Moorcroft  (1963»  and  Benson  (I960). 
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Cause  of  Radio  Star  Fadeout 


It  can  be  shown  (Moorcroft  1963,  Flood  1963)  that  the  fringe  amplitude 
output  of  a  phase -switched  (or  phase- swept  interferometer  of  the  type 
described  by  Lawrence  and  Jespersen  1961)  is  given  by 


B-(l) 


where  is  the  power  emitted  from  the  radio  source 

is  the  mean  squared  amplitude  fluctuation  of  the  star  signal 

is  the  spatial  correlation  function  of  the  amplitude  fluctuations 

is  the  mean  squared  phase  fluctuation  introduced  by  a  random 
ionosphere 

^Cr)  is  the  spatial  correlation  function  of  the  phase  fluctuations 

In  the  derivation  of  the  above  formula,  the  phase  fluctuations  have 
been  assumed  to  be  normally  distributed  with  zero  mean  value  (a  reasonable 
assumption)  and  the  amplitude  and  phase  fluctuations  have  been  assumed  to 
be  independent  (a  not  so  reasonable  assumption). 

Now  the  most  intense  amplitude  scintillation  phenomena,  usually 
associated  with  large  phase  fluctuations,  are  such  that  the  star  signal  ampli¬ 
tudes  at  the  two  antennas  are  Rayleigh-distributed  and  totally  uncorrelated. 
Under  these  circumstances 


==- * &nd 


so  that 


B-(2) 
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The  maximum  reduction  in  interference  fringe  visibility  due  to 
amplitude  scintillations  alone  is  21  percent  or  roughly  1  db.  Since  the 
fadeouts  which  form  the  subject  of  this  report  show  reductions  in  fringe  ampli¬ 
tude  of  2  db  to  greater  than  18  db,  it  is  clear  that  fadeouts  are 
primarily  caused  by  phase  fluctuations.  The  relationship  between  the  depth 
of  the  star  fade  and  the  mean  squared  phase  fluctuations  observed  on  a 
baseline  of  length  r  is  approximately 

db  decrease  =  ^-3  B-(3) 

Further  evidence  that  phase  fluctuations  are  the  prime  cause  of  fade- 
outs  can  be  adduced  from  Figuze  B-8 where  simultaneous  observations  at  52 
and  147  Mcs  of  star  fades  between  2  and  10  db  (the  limit  of  good  accuracy 
is  reading  the  data)  are  plotted.  The  mean  values  of  the  data  are  seen  to 
be  fitted  with  a  A*  law  which  is  just  what  one  would  expect  for  phase 
fluctuations  due  to  electron  density  fluctuations  in  the  ionosphere  (Booker 
1958). 


Anisotropy  of  the  Irregularities  Producing  the  Phase  Fluctuations 

The  analysis  so  far  has  dealt  with  data  from  the  East-West  baselines. 
The  simultaneous  operation  of  North-South  and  East -West  baseline  offers 
the  opportunity  to  examine  the  anisotropy  of  the  correlation  function. 

By  restricting  the  analysis  to  those  fadeouts  which  occur  within  one 
hour  of  the  meridian  passage  of  the  radio  source,  one  can  compare  the  db 
loss  on  the  North-South  baseline  with  the  db  loss  on  the  East- West  baseline. 
In  making  these  comparisons,  the  effective  North-South  baseline  length 
is  the  physical  length  multiplied  by  the  sine  of  the  elevation  angle  (10**  and 
74*  respectively,  for  lower  and  upper  culmination  respectively). 

Not  all  fadeouts  occurring  during  meridian  passage  are  included  in  the 
analysis  which  follows.  These  cases  of  total  star  fadeout  obviously  could 
not  be  analyzed.  Decisions  regarding  the  relative  strengths  of  the  fadeouts 
on  the  two  baselines  were  possible  in  124  cases.  Of  these,  85  cases  show/ed 
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the  fadeout  to  be  atronger  on  the  £ast>Weat  baseline  and  39  cases  indicated 
that  the  fadeout  was  equal  or  stronger  on  the  North>South  baseline.  Almost 
all  of  the  124  cases  analyzed  occurred  during  the  time  of  lower  culmination 
when  the  line  of  sight  from  the  interferometers  to  the  radio  source  is  nearly 
perpendicular  to  the  earth's  magnetic  field.  The  fact  that  39  events  imply 
a  shorter  correlation  distance  along  the  direction  of  the  magnetic  field  is 
most  surprising  since  all  reports  of  star  scintillation  experiments  indicate 
that  the  correlation  distance  (of  the  electron  density  fluctuations  causing 
the  star  scintillations)  measured  along  the  magnetic  field  is  something  like 
five  times  the  correlation  distance  measured  transverse  to  the  magnetic 
field  (Spercer  1955). 

A  partial  explanation  of  the  apparent  discrepancy  between  our  results 
and  previous  reports  of  the  anisotropy  of  the  correlation  function  may  lie 
in  the  fact  that  the  normal  scintillation  measurements  apply  to  F-region 
irregularities.  The  auroral  dependence  previously  indicated  might  imply 
that  some  of  the  fadeouts  here  reported  may  actually  be  caused  by  the 
auroral  E  region  where  the  anisotropy  of  the  large  scale  irregularities  is 
not  well  developed. 

Correlation  of  this  anomalous  anisotropy  of  the  phase  fluctuations 
with  K  figures  from  Agincourt  has  been  attempted.  The  average  K  figure 
for  those  fadeouts  which  showed  a  greater  signal  loss  on  the  East>West 
baseline  was  2.  39 >  while  for  those  fadeouts  which  showed  a  greater  loss  on 
the  North-South  baselines  the  average  K  figure  was  2.  47.  It  would  appear 
that  magnetic  activity  has  no  strong  correlation  with  the  anisotropy  of  the 
phase  fluctuations.  It  should  be  noted  however  that  only  six  instances  of  a 
fadeout  at  upper  culmination  were  considered  in  this  analysis  (all  of  which 
were  stronger  on  the  East-West  baseline)  so  that  the  other  95%  of  the  data 
corresponds  to  lower  culmination  when  the  E  and  F-region  ionospheric 
penetration  points  are  350  and  900  kilometers  north  of  Agincourt,  the  source 
of  the  K  figure  determinations. 
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Correlation  'vith  Spread*F  Activii 


The  F  plots  published  by  the  Defense  Research  Communications 
Establishment  (Canada)  for  the  ionosonde  in  Ottawa,  Ontario  have  been 
examined  for  the  diurnal  variation  of  Spread-F  at  Ottawa  during  the  period 
September  1,  I960  through  August  31,  1961.  All  half  hourly  readings 
indicating  a  spread  were  tabulated  and  the  results  are  plotted  in  Figure  B-9. 
As  expected,  Spread-F  is  a  nightti*  e  phenomenon,  occurring  88  percent  of 
the  nights  at  0300  local  time.  The  location  of  the  maximum  of  Spread-F 
occurs  approximately  six  hours  after  the  evening  maximum  of  star  fades. 

The  fact  that  the  maximum  in  the  diurnal  variation  of  star  scintillations 
precedes  the  Spread-F  maximum  has  been  reported  by  numerous  investigator 
and  the  time  lag  between  the  two  phenomena  :s  sometimes  ascribed  to  the 
shielding  of  higher  altitude  irregularities  which  presumably  cause  star 
scintillations  by  the  lower  ionosphere. 

The  fact  that  two  processes  indicate  a  nighttime  maximum  may 
introduce  a  correlation  into  processes  which  are  otherwise  independent. 

With  only  one  year's  data  on  hand  we  do  not  have  a  sufficient  number  of 
cases  to  look  for  correlation  in  the  fluctuations  of  each  process  from  their 
own  average  diurnal  variation.  It  is  interesting  to  note,  however,  that  the 
300  kilometer  peaet  i'ation  point  at  the  time  of  lower  culmination  is  above 
James  Bay  and  just  south  of  the  permanent  maximum  of  Spread-F  reported 
by  Penndorf  (1962)  for  the  Foxe  Basin. 

Size  of  Irregularities 

The  size  of  the  irregularicies  producing  star  fades  can  be  deduced  from 
experimental  data  and  Equation  B(4.  R  should  be  noted  that  in  order  to  plot 
the  phase  correlation  function,  the  magnitude  of  ,  the  mean  squared 

phase  iluctuations  for  two  totally  uncorrelated  paths,  must  be  known.  The 
experimental  data  on  hand  are  limited  to  a  maximum  signal-to-noise  ratio 
of  approximately  20  db  and  since  a  loss  of  20  db  on  a  900  foot  East-West 
baselines  is  not  uncommon,  ^he  estimates  which  follow  are  only  order  of 
magnitude.  It  would  appear  from  our  analysis  that  a  maximum  scale  size 


PEftCENTAQE  HOURS  SNOViHS  SPREAD 


LOCAL  STANDARD  TINE 


Figure  B-9  PER  :ENT  OF  HOURS  SH0MIN6  SPREAD  F  AT  OTTAWA,  ONTARIO 
SEPTEMBER  I,  I960  <  AUGUST  31,  1961 
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of  200>300  meters  is  not  inconsistent  with  data  taken  during  intense  star 
fades.  This  figure  is  deduced  by  noting  that  the  star  fade  loss  increases 
rapidly  between  the  600  and  900  foot  spacings  implying  that  there  is  a 
rapid  decr^^a^e  in  the  correlatir’*  •'••aluatcd  zt  spacir.^;  cT  CCC  and 

900  feet.  It  is  recognized  that  an  estimate  of  this  sort  probably  tends  to 
underestimate  the  size  of  the  irregularities.  We  have  not  attempted  to 
estimate  the  correlation  distance  in  the  North-South  direction  since  we 
have  <mly  two  baselines  and  therefore  two  data  points.  In  addition,  since 
the  vast  majority  of  the  fadeouts  occur  near  lower  culmination,  our  effective 
baseline  dimensions  in  the  North-South  direction  are  only  190  and  380  feet. 


Radio  Star  Fadeouts  in  Hawaii 


Interferometer  recordings  of  the  radio  sources  in  Cygnus  and  Cassiopeia 
were  obtained  on  the  island  of  Maui  (latitude  21*N)  during  the  period  May  15, 
1962  through  November  5,  1962.  Baselines  of  150,  300  and  600  feet  East-West 
were  used  at  52,  147  and  253  Mcs  while  300  and  600  foot  North-South  baselines 
were  used  on  the  two  lower  frequencies.  During  this  perioa,  one  definite 
radio  star  fadecut  was  recorded,  beginning  at  1020  hours  GMT  on  July  19,  1962 
and  lasting  until  roughly  >800  hours  GMT. 

The  field  site  was  almost  unsuitable  for  scintillation  work.  Quasi- 
continuous  power  line  arcing  created  a  very  serious  ir^terference  problem 
so  that  the  exact  magnitude  of  the  fadeout  could  not  be  measured.  However, 
the  significant  features  of  the  records  are  as  follows: 

1.  The  intensity  of  the  signal  loss  was  greater,  the  lower  the 
operating  frequency  and  the  greater  the  baseline  line  lengths; 

2.  The  duration  of  the  radio  star  fadeout  increased  with  increasing 
baseline  length  and  decreasing  operating  frequency; 

3.  The  duration  and  intensity  of  fadeout  on  the  North-South  baselines 
was  smaller  than  the  duration  and  intensity  on  the  East-West 
baselines. 

4.  The  fadeout  was  complete  -  that  is,  the  signal  dropped  to  the 
noise  level  for  at,  least  ten  minutes  on  the  2^3  Mcs  bOO  foot 
East-West  baselin**. 
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During  this  time  period,  H,  F.  signals  from  Midway  Island  were 
monitored.  The  quality  of  the  signals  was  characterized  by  an  unusual 
flutter  fade.  The  available  data  suggest  that  this  was  an  instance  of  the 
equatorial  F>  region  flutter  fade  reported  by  Koster  (1963).  The  fact  that 
there  was  cmly  one  instance  of  this  fairly  common  equatorial  phenomenon 
may  well  be  explained  by  the  geomagnetic  latitude  of  Maui  (21  *N).  It  is 
also  worth  noting  that  the  data  are  consistent  with  field  aligned  irregularities 
having  a  high  axial  ratio. 

400  Mcs  Scintillation  in  Trinidad 

A  400  Mcs  radiometer  employing  an  eighty-five  foot  diameter  paraboloid 
monitored  the  radio  sources  in  Cygnus,  Cassiopeia  and  Centaur  us  during  the 
night  of  July  15,  1960.  (These  data  were  taken  by  the  General  Electric 
Company  and  were  made  available  by  the  Rome  Air  Development  Center. ) 
Prior  to  2107  hour  EST,  the  signads  from  Cygnus  were  normal,  indicating 
a  0.  1  db  rms  fluctuation.  An  abrupt  start  of  scintillation  activity  was  noticed 
at  2107  hours;  there  is  an  uncertainty  of  the  exact  time  of  the  start  of 
scintillation  activity  due  to  a  minor  equiprnent  malfunction.  Figure  B-10 
example  of  the  scintillation  of  the  source  in  Cygnus.  The  zero  signal  level 
is  at  an  ordinate  value  of  1.  4  and  the  mean  signal  level  during  imdisturbed 
periods  has  an  ordinate  value  of  5.  (The  ordinate  scale  is  linear  in  power. ) 
Figure  B-11  taken  some  time  earlier,  is  a  drift  plot  of  the  Cygnus  source 
through  the  radiometer  beam. 

During  the  course  of  these  observations,  the  antenna  was  also  turned 
to  observe  the  radio  source  in  Centaurus  in  the  Souih-West  and  Cassiopeia 
in  the  North-East.  The  Centaurus  signals  were  normal--  no  pronounced 
scintillation  activity  was  discernible.  The  Cassiopeia  signal  on  the  other 
hand  fluctuated  strongly. 

Since  the  recordings  are  linear  in  power,  scaling  of  the  records  for 
the  fluctuation  of  the  power  AP  from  the  mean  value  of  the  power,  P  , 
was  straightforward.  Thble  Bill  lists  some  of  the  values  of  the  root  mean 
square  of  APfP  which  were  measured. 
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Figure  B>IO  SCINTILLATION  OF  RADIO  SOURCE  IN  CYQNUS,  OBSERVED  AT  400  HCS 
FROM  TRINIDAD  B.W. I.  -  JULY  15,  I960 
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Figure  8-il  DRIFT  PATTERN  OF  RADIO  SOURCE  IN  CYONUS 
JULY  15,  I960  -  IKK)  NCS  -  TRINIDAD 
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Table  B-IH 

RMS  Floctuation  in  Power  15  July  1960 


Time 

Radio  Source 

2204-2213 

Cassiopeia 

0.24 

2239-2249 

Cygnus 

0. 16 

2315-2330 

Cygnus 

0.39 

2400-0010 

Cygnus 

0.52 

The  highest  percentage  fluctuation  measured  was  0.52.  It 
should  be  noted  the  percentage  fluctuation  of  power  for  a  Rayleigh  distribution 
of  amplitudes,  which  would  imply  the  absence  of  an  undiffracted  specular 
wave,  is  ^00%.  Clearly,  then,  under  the  circumstances  of  July  15,  I960, 
multiple  scattering  was  not  complete  and  an  undiffracted  component 
of  the  signal  was  present. 

These  Trinidad  data  are  decidedly  different  from  the  usual  data 
obtained  at  this  site  and  are  not  meant  to  be  ’’typical.  ”  The  planetary 
K  index  for  this  particular  period  was  at  the  8*^  level  and  we  are  not  sure 
whether  the  data  pertain  to  an  aurorally  disturbed  E  or  F  region. 
Nevertheless,  the  data  are  pertinent  in  that  they;  along  with  the  data, 

represent  unusually  high  values  of  phakc  and  amplitude  scintillations 
measured  at  low  latitudes.  There  is  reason  to  believe  that  the  Trinidad 
data  may  refer  to  auroral  irregularities  but  the  Maui  data  most  probably 
refer  to  F-region  irregularities. 

The  Effect  of  Sunspot  Activity  on  Radio  Star  Fadeouts 

In  an  attempt  to  estimate  the  effects  of  solar  activity  upon  the 
occurrence  of  radio  star  fadeouts,  the  52,  147,  and  253  Mcs  interferometers 
were  operated  from  Buffalo,  N.  Y.,  during  the  period  December  1,  1963 
through  May  31,  1964.  The  intensity  and  duration  of  radio  star  fadeouts 
is  apparently  much  less  now  than  during  a  comparable  period  in  1960-61. 

It  would  appear,  therefore,  that  the  number  as  well  as  strength  of  the 
radio  star  fadeouts  is  associated  with  solar  activity. 
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Star  Fadeouts  and  Resolution 


To  obtain  a  semimaximuni  value  for  ,  we  have  to  examine  the 

data  carefully.  There  have  been  a  large  number  of  occasions  when  the 
52  Mc/s  interferometers  showed  a  complete  loss  of  star  signal  on  all  base 
lines.  Since  the  maximum  signal*to*noise  ratio  is  of  the  order  of  20  db, 
an  18-  *.o  20-db  star  fade  is  analogous  to  complete  equipment  saturation. 

All  we  can  say  is  tlu!.t  the  loss  exceeds  18  db  and  therefore  that  the  mean 
squared  phase  fluctuations  are  in  excess  of  4.  2  square  radians  at  52  Mc/s. 
However,  on  at  least  three  occasions,  a  147  Mc/s  star  fadesut  of  greater 
than  11  db  has  been  recorded.  This  would  imply  an  88-db  loss  at  52  Mc/s 
on  these  occasions. 

Radio  star  fadeout  observations  in  the  auroral  zone  have  been  made 
by  Little  (1959)  who  reported  that  a  7-db  loss  was  measured  at  a  frequency 
of  223  Mc/s  on  an  interferometer  baseline  of  100  meters.  Using  the 
law,  we  would  predict  a  loss  of  16  db  at  147  Mc/s  rather  than  the  11  db  used 
in  the  preceding  paragraph.  Since  the  1 1  db  was  not  an  absolute  maximum 
and  we  would  expect  stronger  auroral  disturbances  at  the  Alaskam  location 
of  Little,  his  observations  are  in  accord  with  the  basic  auroral  nature  of  the 
phenomenon  as  well  as  the  validity  of  the  /f  dependence.  From  Equation  B-3 

l^(r)  —  (qp— iy)db  less 

^  30  square  radians  m,i  SO  ntcs  B-(4) 

(using  Little’s  data).  The  mean  squared  phase  fluctuations  at  any  other 
frequency  will  therefore  be  assumed  to  have  a  "maximum"  (in  view  of  the 
A*  dependence). 


square  radians 
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As  typical  examples  of  the  variance  of  the  phase  fluctuations,  we  would 
expect  36  x  10’^  and  9  x  10"'*  square  radian  at  frequencies  of  5000  and 
10,  000  Mcs  respectively..  That  these  values  will  net  be  uncommon  can  be 
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All  we  can  say  is  tluit  the  loss  exceeds  18  db  and  therefore  that  the  mean 

I 

squared  phase  fluctuations  are  in  excess  of  4.  2  square  radians  at  52  Mc/s. 
t  However,  on  at  least  three  occasions,  a  147  Mc/s  star  fadesut  of  greater 
'  than  11  db  has  been  recorded.  This  would  imply  an  88-db  loss  at  52  Mc/s 
on  these  occasions. 

Radio  star  fadeout  observations  in  the  auroral  zone  have  been  made 
by  Little  (1959)  who  reported  that  a  7-db  loss  was  measured  at  a  freqviency 
of  223  Mc/s  on  an  interferometer  baseline  of  100  meters.  Using  the  Af 
law,  we  would  predict  a  loss  of  16  db  at  147  Mc/s  rather  than  the  11  db  used 
in  the  preceding  paragraph.  Since  the  11  db  was  not  an  absolute  maximum 
and  we  would  expect  stronger  auroral  disturbances  at  the  Alaskan  location 
of  Little,  his  observations  are  in  accord  with  the  basic  auroral  nature  of  the 
phenomenon  as  well  as  the  validity  of  the  Af  dependence.  From  Equation  B-3 

(^^)db  less 

30  square  radians  «/  SO  ?>tcs  B-(4) 

(using  Little's  data).  The  mean  squared  phase  flucLuations  at  any  other 
frequency  will  therefore  be  assumed  to  have  a  "maximum"  (in  view  of  the 
A*  dependence). 


B-(5) 


As  typical  examples  of  the  variance  of  the  phase  fluctuations,  we  would 
expect  56  x  lO’"^  and  9  x  10"'^  square  radian  at  frequencies  of  5000  and 
10,  000  Mcs  respcctivelv.;  That  these  values  will  net  be  uncommon  can  be 
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ascertained  from  the  fact  that  estimates  of  a  "typical"  strong  fade  ai  ^ 

147  Mc/s  are  approximately  7  db.  This  would  reduce  the  above  "maximum" 
values  by  a  factor  less  than  Z. 

The  above  estimates  of  the  variance  of  the  phase  fluctuations  at  the  ^ 

S  and  X  bands  are  surprisingly  large  smd  indicate  that  for  some  long  ^ 

propagation  paths  which  traverse  the  auroral  ionosphere,  the  ionospheric 
contributions  to  the  short  term  variance  of  the  phase  fluctuations  are  of  the 
same  order  as  the  tropospheric  contributions. 


Effect  of  Ionospheric  Phase  Fluctuations  on  Angular  Resolution 

The  lack  of  phase  correlation  at  two  points  implies  a  loss  of  antenna 
gain  and  angular  resolution.  If  we  define  the  maximum  antenna  dimension  D 
as  that  dimension  wluch  would  result  in  a  mean  squared  phase  fluctuation  of 
1  radian  squared,  then  estimates  of  the  maximum  resolving  power  (based  on 
a  Rayleigh  criterion  of  A/D)  cam  be  made  from  the  experimental  data. 

Following  Ratcliffe's  (1956)  usag**.  we  will  take  greater  tham  one  squ 

radian  to  imply  multiple  scatter,  indie uting  that  the  intensity  of  the  undiffracted 
wave  is  not  much  larger  than  the  intensity  of  the  waves  arriving  over  a  variety 
of  angles.  For  the  case  small  compared  to  1  radian  squared,  the 

scattered  waves  are  small  compared  to  the  zero-order  wave,  and  the  angle 
of  arrival  is  "perfectly  "defined. 

FromB(5)  we  conduce  that  for  frequencies  in  excess  of  300  Mc/s 
the  angular  resolution  will  be  determined  by  the  size  of  the  antenna  aperture. 

At  frequencies  less  than  300  Mc/s,  the.  limits  to  angular  resolution  vail  be 
.l.:termined  by  the  ionospheric  phase  path  irregularities.  The  cone  angle  of 
^  -rival  for  frequencies  less  than  300  Mc/s  is  determined  by  where  D 

is  as  defined  above. 

If  the  phase  fluctuations  are  df  lined  by 

^r)  -  S'  B.(6) 


B-{9) 


squ 

ted 

ty 


Equation  B(9)  is  analogous  to  the  result  of  Ratcliffe  (1956)  who  points 
out  that,  in  diffraction  from  an  irregular  screen,  when  the  mean  squared 
phase  fluctuations  exceed  1  radian  sqtiared,  the  effective  correlation 

distance  measured  on  the  ground  is  equal  to  the  true  correlation  distance 
divided  by  . 

The  cone  angle  of  arrival  6  is  defined  by 
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It  is  clear  that  to  predict  the  amount  of  beam  broadening  to  be  expected 
at  frequencies  less  ihan  300  Mc/s,  a  knowledge  of  the  correlation  distance  L 
is  required.  For  ionospheric  phase  path  fluctuations  at  least,  L  is 
generally  greater  than  100  meters,  since  the  mean  squared  phase  fluctuations 
generally  increased  at  distances  greater  than  300  feet  at  147  Mc/s.  To  put 
a  pessimistic  limit  to  the  degree  of  resolution  that  can  be  achieved,  we 
take  L  to  he  100  meters  and  to  be  equal  to  square  radians 

(equation  B(5) )  Tkble  B-IV  lists  the  limitations  to  angular  resolution  as  a  function 
of  frequency. 
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Table  B-IV 

Limitations  to  Angular  Resolution  as  a  Function  of  Frequency 

Frequency,  Mc/s 

50  150  300 

0  In  milliradians  360  40  1  0 


Bandwidth  Limitations  as  a  Consequence  of  Phase  Path  Fluctuations 

The  observed  relationship  between  radio  star  fadeouts  and  auroial  back- 
scattering  at  closely  allied  frequencies  suggests  that,  on  occasion  at  least, 
the  source  of  the  fluctuations  is  in  the  auroral  E  region.  Under  these  circum¬ 
stances,  the  geometry  in  Figure  B-i2.  is  appropriate  for  estimating  the  multipath 
bandwidth  limitations  imposed  by  propagation  through  the  auroral  ionosphere. 
Consider  the  worst  possible  case:  a  source  at  an  infinite  distan  e  from  the 
scattering  screen.  Radiation  received  on  the  ground  over  a  cone  angle  G 
arrives  with  a  time  delay  (referred  to  the  ray  arriving  in  the  line-of- sight 
direction)  of  T 


Equation B-(ll)  implies  that  a  very  narrow  pulse  would  be  broadened 

into  a  pulse  length  V  ,,  Such  pulse  broadening  would  limit  the  range 

CT' 

resolution  of  a  pulse  radar  to  an  amount  , 

cr  ^  \  ^ 

z  "  z  \  ec  )  '  Jk 

and  from  B-(IO) 

CT  ^ 

~Z  1i>L^  B-(IZ) 
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ELCmiON  VIEW  OF  MIMRAL  SEOMETRr 


0*  *  (I  ♦£?) 

4 


Figure  B- 1 2  AURORAL  MULTIPATH  GEOMETRY 


Takixig  ^  to  be  of  the  order  of  8  x  10^  roetere  and  I,  to  be  100  meters 
(minimum),  we  obtain  a  pessimistic  limitation  to  range  resolution: 


cr 
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Range  resolution  of  better  than  1  foot  (approximately  0.  3  meter)  can  be 
achieved  at  ali  wavelengths  less  than  0.  5  meter  if  multipath  al<me  is  the 
limiting  factor  in  range  resolution. 
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APPENDIX  C 

TROPOSPHERIC  PHASE  PATH  FLUCTUATION  EXPERIMENT 


Late  in  1962  it  became  evident  that  shifts  in  Air  Force  planning 
would  preclude  the  availability  of  a  suitable  satellite  vehicle  for  the 
WIDE  BAND  experiments.  Therafter,  a  short-term  contract  (AF  30(602)- 
2840)  was  opened  to  permit  preliminary  examination  of  the  question:  what 
portion  of  the  satellite -ground  propagation  patf  behavior  can  be  ascribed 
to  tropospheric  influences?  Under  this  contract,  a  B-57  airplane  was 
fitted  with  an  S-band  beacon  transmitter  and  there  were  begun  a  series 
of  flights  during  which  transmissions  were  made  over  an  approximately 
SO-mile  long,  constantly  varying  tropospheric  path  from  air  to  the  inter¬ 
ferometer  ground  stations.  The  aircraft  speed  was  typically  400  knots 
and  its  altitude  42,000  feet. 

A  report  summarizing  the  results  from  the  tropospheric  propagation 
experiment  was  issued  at  the  termination  of  contract  AF  30(602) -2840, 
and  the  experimental  effort  was  resumed  under  contract  AF  30(602)-2077. 
Specific  goals  for  the  continued  program  were  the  investigation  of 
tropospheric  phase  path  fluctuations  under  a  variety  of  different  meteor¬ 
ological  conditions,  and  in  particular  during  thunderstorms.  Because 
of  the  large  quantity  of  data  collected  during  the  Spring  of  1963,  the 
additional  data  to  be  collected  during  the  continuation  of  the  program  were 
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intended  to  represent  the  many  different  weather  conditions  encountered 
in  the  Buffalo  area  rather  than  just  more  data  for  the  average  good  weather 
day.  Because  of  weather  ccmditions  prevailing  during  July  and  August 
1963,  only  two  fli|^ts  were  conducted  during  disturbed  (thunderstorm) 
weather.  The  remaining  data  represent  meteorologically  "quiet" 
conditions. 

EXPERIMENTAL  PROCEDURE 

The  tropospheric  phase  path  fluctuation  experiment  was  conducted 
at  the  CAL  Sand  Hill  receiving  site^^^  located  approximately  20  miles 
northeast  of  Buffalo,  New  York.  The  purpose  of  the  program  was  to 
obtain  experimental  data  that  could  be  used  to  determine  the  limitations 
on  angular  radar  resolution  at  S>band  imposed  by  random  fluctuations  of 
refractive  index  in  the  troposphere. 

In  order  to  obtain  abng  propagation  path  through  the  troposphere, 
the  S-band  (2250  me)  beacon  transmitter  was  carried  on  board  a  B-57 
aircraft  in  the  altitude  range  of  39,000  ft  to  42,000  ft.  The  flight  plan 
followed  by  the  B-57  is  shown  in  Figure  C-1.  This  pattern  was  flown 
twice  during  each  experiment,  and  required  approximately  one  hour  per 
lap  at  an  average  ground  speed  of  400  knots. 


The  beacon  transmitter  carried  by  the  B-57  consisted  of  a 
cavity  oscillator  capable  of  delivering  2  watts  of  c,w.  power  to  its 
antenna.  The  transmitter  and  associated  solid-state  power  supply 
circuitry  were  thermally  stabilized  by  means  of  a  temperature  controlled 
heat  sink.  Using  a  Freon-evaporation  cooling  system  plus  shock  mounting 
of  all  transmitter  components,  resulted  in  long-term  (2  hours)  frequency 
stability  of  about  10  parts  per  million  and  an  incidental  frequency 
modulation  spectrum  that  was  about  4  kilocycles  wide  at  the  halfpower 
points.  The  transmitting  antenna,  a  quarter  wave  stub  against  a 
ground  plane,  was  mounted  on  the  under  side  of  the  fuselage,  near  the 
tail,  with  the  stub  extending  downward.  Occasional  tracking  difficulties 
were  experienced  when  the  antenna  was  shielded  by  the  aircraft  body 
during  a  heavily  banked  turn.  These  did  not  result  in  serious  loss  of 
data. 


The  receiver-interferometers  were  arranged  on  two  lengths  of 
orthogonal  baselines,  lying  north- south  and  east- west,  and  of  lengths 
100  and  1000  meters.  Because  of  the  very  large  (relative  to  the 
operating  wavelength)  baseline  lengths  used  in  this  experiment,  all  local 
oscillators  in  the  interferometer  complex  were  phase  locked  to  a  common 
reference  signal  so  that  phase  fluctuations  due  to  mechanical  expansion 
or  contraction  of  the  interconnecting  cables  were  removed.  The 
V  stalls  of  the  phase  locking  arrangement  used  are  documented  in  a 
report^^^  entitled,  "WIDE  BAND  Receiving  Installation.  " 
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The  analog  data  from  the  interferometera,  proportional  to  the 
cosine  of  the  R.F.  phase  difference  between  the  interferometer  antennas, 
were  recorded  on  magnetic  tape  at  30  in/sec.  A  sample  of  the  analog 
interferometer  output,  reproduced  from  the  magnetic  data  tape,  is 
shown  in  Figure  C-2. 

Whenever  the  aircraft  flew  a  course  parallel  to  one  of  the  1000-meter 
baselines,  very  high  differential  phase  rates  were  observed,  occasionally 
reaching  a  value  of  190  radians  per  second.  In  order  to  obtain  a  least 
count  of  2^  in  phase,  we  digitized  the  analog  data  at  a  rate  of  5000  words 
per  second  into  six-bit  binary  words.  The  digital  data  points  were 
stored  on  magnetic  tape  in  a  format  compatible  with  the  requirements 
of  the  CAL  IBM  704  computer.  This  oomputer  was  used  to  perform  ail 
of  the  subsequent  data  processing. 

In  addition  to  the  four  interferometers,  the  receiving  site  contained 
a  tracking  receiver,  the  purpose  of  which  was  to  follow  the  S-band 
beacon  transmitter  and  transmit  azimuth  and  elevation  information  to 
the  antennas  in  the  interferometer  complex.  The  tracker  antenna  was 
a  16  foot  paraboloid  with  a  conical -scanning  feed.  Tracking  accuracy 
was  10  milliradians  with  wind  speeds  up  to  20  mph.  Continuous  tracking 
information  (sine  and  cosine  of  the  azimuth  angle  and  the  sine  of  the 
elevation  angle),  and  WWV-time,  were  recorded  on  a  four-channel 
Sanborn  recorder  during  each  experiment. 
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Figure  C-2  INTERFEROMETER  SAMPLE  DATA 
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DATA  ANALYSIS  AND  RESULTS 


Theoretical  Approach 

The  purpose  of  the  CAL  experiment  was  to  investigate  the 
magnitude  and  the  spectrum  of  S-band  phase  fluctuations  caused  >y 
random  refractive  index  fluctuations  in  the  troposphere.  Because 
our  experiment  was  conducted  with  the  transmitter  on  board  an  aircraft, 
the  phase  path  between  transmitter  and  receiver  s  was  constantly  changing, 
a  fact  which  made  data  analysis  extremely  difficult.  The  major  problem 
encountered  was  separating  the  phase  changes  due  to  aircraft  motion 
from  those  produced  by  tropospheric  phase  path  fluctuations.  In  the 
discussion  that  follows  the  time  dependence  of  the  phase  changes  due 
to  the  unperturbed  (i.e.,  constant  velocity,  no  vibration,  etc.)  aircraft 
motion  will  be  derived  in  terms  of  aircraft  velocity,  range,  and  the 
angle  between  the  direction  vector  to  the  aircraft  and  the  interferometer 
baseline.  The  difference  between  such  a  deterministic  expression  for 
the  aircraft  motion  and  the  raw  phase  data  must  constitute  the  desired 
random  phase  data  due  to  propagation  through  the  troposphere. 

Figure  C-3  represents  the  static  geometry  of  our  experiment. 

The  transmitter  is  located  at  T,  a  distance  R  from  the  interferometer. 
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Figure  C*3  >  Geometry  representing  transmitter  and  receivers 


The  interferometer  antennas  are  located  at  points  A  and  B  '.m  a  baseline 
that  is  (d)  long  and  oriented  along  the  y-axis.  (The  y»axis  coincides 
with  the  north-south  baseline  of  our  experiment.)  The  data  collected 
d  uring  our  experiment  was  a  function  of  the  phase  difference 
between  antennas  A  and  B,  where 

and  .  The  quantities  and  are  the 

random  phasesincurred  because  of  random  variations  in  the  refractive 
index  along  the  paths  7^-^  and  7^-3  . 
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We  will  next  examine  the  difference  in  magnitode  A  P  between 
vectors  iPn  and  ^  ^  . 


le^l  =f>pv-^ 

Kel*-  e.Z) 


(a) 


(b) 
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The  solution  to^^  from  equation  (C*2)  involves  expc-sding  the  quantities 
l^nl  •”<>  »“  a  power  aeries.  This  results  in  the  following  expression 

for  : 


(e-3) 


AR'  iftil-  /ip6!'/i?l[a--2^%£*4  3£iV -  7 

^  8  It, 

where  Ct  I /  *  Note  that  in  equation  (C- 3) 

ali  the  even  powers  of  j  are  missing.  Thus  if  j  is  a  small  compared 

to  unity,  the  power  series  converges  rapidly  to  seros  In  our  experiment 

the  baseline  lengths  were  equal  to,  or  smaller  than,  1000  meters,  and  the 

range  to  the  aircraft  was  never  less  than  50  km.  Consequently,  by 

neglecting  all  terms  containing  I4l  to  powers  greater  than  unity,  our 

error  will  be  less  than  1  part  of  10,000  in  the  determination  of  R 

an  error  less  than  5  milliradians  when  ^  R  is  multiplied  by 
-  — 

(  is  equal  to  15  //  radians/meter  in  oui  experiment.)  A. phase 

A 

error  of  this  magnitude  is  estimated  to  be  at  least  a  factor  of  ten  smaller 

e 


than  the  rms  val  le  of  equation  (C>1)  and  can,  therefore. 


be  neglected. 
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To  accoont  for  oar  moriaf  aircraft  we  moat  examine  ftio  aame 
geometry,  but  let  the  range  vector  be  time  dependent  ao  that 


Neglecting  parallax,  the  range  difference  become  a 

AROi)-^  \£«^IaU) 


where 


-*  .  A 


(e.-i 

l£«)l  ^ 


(b) 
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The  aubacript  'i  0  )  denotea  the  value  of  R  at  time  \  /Mt  '•  *be 
difference  between  time  zf  o  and  time  • 

It  should  be  evident,  that  the  solution  of  will  once  more 

involve  a  power  eeriea,  this  time  in  powera  of  Thua,  we  obtain 
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VAt 

J^o 


j  S(Po*v)  •  /OfPt  •y)\  ‘  •  J 


The  expansion  can  be  carried  out  to  powers  ci  higher  than 

the  cubic;  the  choice  of  how  far  to  carry  this  process  will  depend  on 
the  magnitude  of  J^"6^  which  is  the  time  over  which  the  aircraft 
can  reasonably  be  expected  to  fly  a  straight  and  unperturbed  course. 
In  general,  we  will  neglect  parallax  and  write  as  a  power 

series  in  time  with  coefficients  so  that 

W 

h»o 


where,  for  the  approximations  made  here,  the  first  four  coefficients 


are: 

a^ldl(h^) 


-/dllljpo-iJXe  •  v)-  (v.^)J 


(b) 


X 


(C) 
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(C-S) 
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It  slMkl  hm  b0I«4  IImM  tfds  power  •orios  ospoaoHm  for  ^  ^ 
cootafair  dot  ooit  vector  Q  wkldi,  io  Figere  C-3,  coiocMcs  with 
the  y-ane.  Thea,  ia  dde  aaalyeie  dM  power  seriee  coeiBcieata  O 
are  correct  for  the  aorth-eoedi  aad  eaat-weet  baaellaee  whea  the 


vector#  ^  ^ad  ^6  j  reepectively.  are  eebetltated  for  ^  . 


hlachtae  Compotattoae 

Aaalof  data  ofacalaed  from  die  experimeat  were  reco.  ded 
oa  magaetic  tape  aad  sabaequeatly  qeaatiaed  iato  six-hit  hiaary  words 
Even  though  the  fuactioa  of  iaterest  to  os  is  die  raw  data 

collected  at  oar  field  site  were  of  die  form  (See  Figure  C-21t 

F0f)‘ 

FM  ^  K^os  I  fM  ^ 


where  ^  (i)  s  O,  i  i,  t  -t  Ct  ^ -t - 

and  are  the  random  phase  data  that  we  wish  *  ->  isolate  and 

examine.  Designate  the  times  at  which  ••  i  • 


(C-9) 


Then 


C-i2 


(b)  (C-iO) 

- . 

(c) 


where  s  2L 


U,  in  equation  (C-10c)«  the  aircraft  position  were  known 
exactly  at  time  4  »  we  could  evaluate  the  polynomial  for 

and  thus  solve  for  A^^^-^^j^the  desired  data  from  this  experiment. 

Because  the  airplane's  position  could  not  be  determined  to  the  accuracy 

s 

commensurate  with  other  data,  we  estimated  values  for  the  coefficients 
in  the  polynomial  for  J  by  computing  a  least>squares  fit  for  the 
raw  phase  data.  In  so  doing,  we  assumed  that,  over  a  five-second  time 
interval,  a  polynomial  of  the  form  'f  would  adequately 

account  for  the  phase  changes  due  to  aircraft  motion.  It  must  be  noted 
that  the  coefficients  R  j  S ,  ^nd  will  differ  from  the  coefficients 

I  I 

Ct  f  L  j  C  in  equation  (C- iOc)  because  the  least-squares 


See  Appendix  C-2  for  a  derivation  of  the  exactness  requirements. 
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p;  oeeBB  docs  not  differcntistc  between  coatribotions  to  the  total  phase 
4oe  to  ttic  aircraft  and  doe  to  the  atmosphere.  Within  tfic 
limitations  of  dhe  assomptions  and  approximations  made,  the  difference 
between  die  total  phase  (€i)  tatd  the  polynomial  f 

is  a  phase  diat  is  caosed  by  random  variations  of  tbs  refractive 
index  along  dM  propagation  path. 


Specific  computations  performed  were  as  follows: 

1.  Compute  die  average  time  between  consecutive  sero  crossings 
of  •  Thi.».(  erage  half  period  is  and  is  defined  as 

T- 

Ca'/J 

where  is  the  time  of  die  first  sero  crossing,  and  6^  ir  die  time 

of  die  last  sero  crossing  in  the  data  interval  of  5  seconds. 

Z.  Compute  the  ilifference  ^  between  the  time  interval 
and  V  •  Normalise  this  difference  to  the  average  half  period  7X 

Then 


(C-IZ) 

Pk-  (&->)-  (  (b) 
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3.  The  values  of  were  next  fitted  with  the  polynomial 

by  the  method  of  least  squares.  The  difference  between 
the  polynomial  (evaluated  at  time  )  and  was  the  remainder 

which  constitutes  the  output  from  the  computer  program.  The  relation¬ 
ship  between  the  polynomials^  ^  ^ "^d  ^ ^  (/iiGii-fC/e'') 

as  well  as  the  similarity  of  developed 

in  the  following  identities: 


Substitute  (C-12b)  for 


>  (b)  (C-i3) 


multiply  both  sides  by  // 


//ei--  ffT-  fl- B€i,- 


Diacussion  of  Results 

On  the  basis  of  the  foregoing  discussion,  and  on  the  assumption 
that  phase  changes  due  to  aircraft  motion  over  a  time  interval  of  five 
seconds  can  adequately  be  represented  by  a  quadratic,  one  would  expect 
that  the  coefficients  Q  )  ^  ^  and  O'  in  equation  (C-lOc)  are  fairly  well 
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approximated  by  the  coefficients  ^  and  0  in  (Cl 3c).*  Recognising 

the  limitations  inherent  to  our  analysis  technique,  we  computed  values 
77^^^  from  six  consecutive  five-second  long  data  sections,  i.e., 
every  five  seconds  new  coefficients  of  fit  (  |^  ,  ^  j  and  C, )  were  computed. 
Data  collected  on  four  different  days  were  analysed.  Of  these,  two  days 
represented  meteorologically  "quiet"  days,  while  shower  and  thunder¬ 
storm  activity  prevailed  in  the  Buffalo  area  on  the  other  two,**  Power 
spectra  were  computed  from  the  values.  These  are  presented 

in  Figures  C-4,  C>S,  C*6,  and  C-7.  Each  of  these  four  figures  is  a 
composite  of  six  40  point  spectra;  each  spectrum  was  computed  from 
5  seconds  of  interferometer  data  from  which  the  polynomial 
had  been  removed. 


The  only  data  available  for  comparison  with  our  results  were 
obtained  with  fixed  propagation  paths.  Under  the  assumption  that 
spectral  densities  obtained  on  fixed  paths  will  scale  directly  with  the 
length  of  the  propagation  path  and  with  the  square  of  the  frequency,  the 
spectral  in  Figures  C-4,  C-5,  C-6  and  C-7  are  in  agreement  (order  of 
magnitude)  with  the  L-band  results  presented  by  Herbstreit  and  Thompson 
in  1955.  In  comparing  their  spectrum  with  ours,  we  have  made  the 
additional  assumption  that  the  frequency  dependence  of  JL  continues 


(1) 


See  Appendix  C-3 

See  Appendix  C-4  for  weather  synopses  and  aircraft  positions  at 
the  times  represented  by  the  data. 


V/ 
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S(f)  RAO^/cm 


Figure  C-l  PHASE  FLUCTUATION  SPECTRUM  I  AUGUST  1963,  1507  EOT 

1000  METER  NORTH- SOUTH  BASELINE 
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*(f)  tiO*/cM 
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Figure  C-7  PHASE  FLUCTUATION  SPECTRUM  28  AUGUST  1963.  1623  EOT 

1000  METER  NORTH-SOUTH  BASELINE 


beyond  the  data  points  presented  by  them.  The  validity  of  this  assumption 
is  substantiated  by  the  9,400  me  phase -variation  spectrum  published  by 
Thompson  and  Janes  in  1959. 

There  are  several  features  common  to  all  of  the  curves  presented 
here.  Probably  the  most  striking  one  is  the  occurrence  of  large,  isolated 
spectral  densities  at  frequencies  near  10  cps,  15  cps,  and  20  cps.  We  do  not 
believe  that  these  spikes  are  caused  by  propagation  through  the  atmosphere; 
the  most  likely  explanation  for  their  existence  is  that  of  wind-induced,  natural 
resonant  vibration  of  either  the  antenna  towers,  dishes,  feed  supports  or  a 
combination  of  all  of  these.  These  spikes,  therefore,  should  be  thought  of 
as  defining  the  system  noise  level. 
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A  second  feature  common  to  practically  all  of  the  spectra  obtained 
from  our  data  is  their  similarity  to  the  spectrum  of  the  time  function  g(t)  =  at^.* 
This  similarity  is  particularly  noticeable  on  meteorologically  "quiet"  days 
(See  spectra  for  August  22  and  27).  Why  this  parallelism  should  exist  and 
what  significance  should  be  attached  to  it.  is  not  well  understood.  It  should 
be  pointed  out,  however,  that  for  power  spectra  to  be  alike,  the  only  require- 
ment  is  that  their  autoco'-relation  functions  be  alike.  This  statement  implies 
nothing  concerning  the  similarity  of  the  time  functions  giving  rise  to  these 
spectra.  To  illustrate  this  point,  we  have  plotted  in  Figure  C-8  the  time 
history  of  (equation  (Cl 3c))  for  a  5 -second  time  interval  on 

28  August  1963.  from  16*'23”'38  ***■  to  16^23”'43***^.  In  the  same  figure 
appears  a  dotted  curve  which  represents  the  time  function  at^  from  which 
the  mean  value  has  been  removed,  and  which,  over  a  5-second  time  interval 
has  the  same  variance  as  the  time  function  .  There  is  no  readily 

apparent  similarity  between  these  two  time  functions. 

In  general,  our  spectra  over  the  frequency  range  from  0.6  cps  to 
20  cps  are  in  qualitative  agreement  with  extrapolated  spectra  obtained 
from  fixed  propagation  paths  and  different  operating  frequency 
However,  we  have  shown  in  Appendix  C-3  that  the  cubic  term  (in  the  poly¬ 
nomial  which  accounts  for  phase  changes  due  to  aircraft  position),  also 
gives  rise  to  a  spectrum  having  the  expected  spectral  densities  and  the 
expected  frequency  dependence.  Therefore,  our  results  must  be  tieated  as 
inconclusive  as  far  as  dete  .'mination  of  the  true  magnitude  and  spectrum  of 
the  phase  path  fluctuations  is  concerned. 


See  also  Appendix  C-3 
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POLINOMiNAL  HkS  BEEN  REMOVED  AS  A  FUNCTSON  OF  TIME 


CONCLUSIONS  AND  RECOMMENDATIONS 


The  tropoapheric  experiment  conducted  by  thia  Laboratory  waa 

intended  to  yield  data  from  which  the  mean  aquare  phaae  error  and  the 

power  apectrum  of  phaae  fluctuationa  due  to  tropoapheric  inhomogeneitiea 

could  be  computed.  It  waa  expected  that  the  mean-aquered  phaae  error 

cauaed  by  S-band  propagation  through  turbulent  media  (such  as  thunder- 

atorma)  would  be  appreciably  larger  than  the  mean-aquare  phase  errors 

arising  from  a  meteorologically  quiet  medium.  This  expectation  is  not 

supported  by  our  data.  Wr  find,  however,  that  the  total  mean- square  phase 

error  due  to  the  combination  of  atmosphere,  aircraft  motion,  antenna 

vibration  and  thermal  noise  is  never  larger  than  10  ^  square  radians 

over  a  time  interval  of  five  seconds.  We  can  safely  conclude,  therefore, 

that  the  phase  error  due  to  atmospheric  inhomogeneities  alone  is  smaller 

than  this  value.  One  of  the  goals  of  the  experiment  was  to  ascertain 

the  limitations  imposed  by  the  troposphere  on  angular  resolution  at 

S-band.  We  can  positively  state  ttat  at  all  times  angular  resolution 

-5 

of  better  than  6  x  10  radians  is  achievable  while  propagating  through 
the  troposphere.  This  estimate  constitues  a  pessimistic  bound,  and 
actual  resolution  limits  should  be  finer. 

Our  experimental  spectra  are  in  qualitative  agreement  with  extra¬ 
polated  results  obtained  by  the  National  Bureau  of  Standards.  Because 
these  spectra  also  have  the  same  form  as  would  be  expected  from  residual 
aircraft  motion  effects,  i.e.,  the  at^  term,  we  must  treat  our  results 
as  being  inconclusive. 
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Because  the  achievement  of  niaximum  system  accuracy  necessitates 
that  angular  value  of  line  of  sight  be  known  to  lO^radiana,  there  should  be 
no  further  pursuit  of  such  propagation  problems  by  airborne  experiments  in 
which  that  precision  cannot  be  realized.  This  recommendation  is  pessismistic 
in  that  it  relates  to  phase  determination  within  0.5  degree.  Relaxation  of 
that  limit  would  be  accompanied  by  relaxation  of  system  precision  requirements. 
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Appendix  C-1 


Detail*  of  Data  Analyeie  Program  Initially  Used 
in  an  Effort  to  Separate  Vehicle  Motion  From  the 
Atmospheric  Effects  on  the  Phase  of  a  Radio  Wave 

During  the  early  stages  of  the  data  collection  and  anal  sis  program 
it  was  recognised  that,  in  order  to  accurately  assess  the  effects  that 
the  atmosphere  has  on  the  phase  of  a  radio  wave,  the  distance  between 
the  radio  source  and  the  receiving  instrumentation  would  have  to 
be  extremely  well  known.  In  fact,  the  precision  to  which  one  has  to 
know  this  distance  is  such  as  to  be  almost  unattainable  even  by  the 
most  modern  methods.  In  an  attempt  to  circumvent  these  extreme 
requirements,  we  stipulated  that  over  the  time  interval  required  for 
the  phase  path  to  change  by  one  wavelength,  the  aircraft  must  follow 
a  straight-line  path  for  which  perturbations  in  motion  would  cause 
phase  fluctuation  errors  that  are  negligible  compared  to  the  propagation 
effects.  On  the  basis  of  this  assumption  alone,  we  proceeded  to  analyse 
the  phase  fluctuations  data  along  the  following  lines  of  reasoning. 

The  phase  difference  between  two  antennas  separated  by  a 
distance  (U' ),  both  receiving  a  signal  from  a  transmitter  a  distance 
away,  will  be  (see  Figure  C-1-1) 


-1(1) 
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where 


and  where  (j  OtA-d  ^r  Q  random  phase  fluctuations 

incurred  along  paths  and  because  of  random  fluctuaticms  in 
the  refractive  index  c-long  these  paths. 


Figure  C-1-1 
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In  the  interferometer  experiment  performed  the  point  r 
wae  in  conaUnt  motion,  and  the  analog  output  from  the  interferometer 
waa  of  the  form 

r(  O-  k  j^f 

where  f  and  »»<> 

where  it  has  been  assumed  that  The  purpose  of  the 

data  analysis  program  is  separate  the  effect  of  the  aircraft  motion 
represented  by  )  from  the  random  phase  • 


(a) 

C-l(2) 

(b> 
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If  tlie  value  of  ('f)  is  small,  then  equation  |C>i  2b)  can  be 
escpanded  to  yield 


(f- f{^)  ^  <f(f)  fcoj  w 


which,  for 


becomes 


(r.4(3) 


■  FC^)  ^  KsU^  f(f)  /  k:  fM 


One  can  reason,  therefore,  that  when  multiple 

of  then  KJ  is  at  its  peak  value.  This  occurs  at  times  *  *  ' 

If  one  can  next  assume  that  f)  varies  linearily  with  time  over  a 
half  period,  then  at  time  ^2  •“  ^  ^  ■^°'*^**  equal  a  multiple 


of  //  radians,  and  Consequently,  the  desired  ’ 

phase  fluctuation  data  should  be  available  when  the  functional  value  of 
is  evaluated  at  a  time  that  is  halfway  between  the  occurrence  of  the 
peak  values  of  )  • 
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We  will  next  re-examine  thia  data  analysie  procedure  and 


show  that  the  conclueiona  reached  arc  not  correct. 

Again,  let  .  At  time  ti 

(Figure  C-1-2), 

F  ad  *  F'  (a) 

{U>)  *  --r 

(b) 

At  time 

FC^s)  --  -/<  (») 

(b) 


C-l(4> 


C-l(5) 


C-l(6) 
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Now  assume  that  varies  linearity  with  time.  Then 


f ^<4/3 -jr- 1 4 (€.) / J 


c-i{7r 


which  for 


become. 


4(^) 


cU 


^  L 


C-l(8) 


Let  the  half  period  be  7^  -  t^3  - 

V 


.  Then 


4l  i  61  . 


-  ^  // 


0-1(9) 
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Let  the  component  of  the  frfiaee  mt  angular  frequency  (69)  be  given  by 
exp  J  Go  t  =  rewrite  equation  C-l(8b)  as 


F«i  3  -  -  ->>tr  Y  ^ j  T'^) 

r({i ) '  Y  y)  ) 


(b) 


c 


(c> 


L 


(d) 


Equation  C-l(lOd)  is  the  actual  output  from  the  data  analysis 
program.  It  differs  fA'om  the  desired  result  in  that  it  contains  a  filter 
function  proportional  to  ^  ^  ^  "Y  *  ^bere  JL^  is  the 

frequency  of  the  analog  data  (interferometer  output).  This  filter  function 
severely  attenuates  the  low  frequency  components  contained  in  • 

Its  crossover  frequency  occurs  at  ^  =  0.805  and  full  traasmission 

occurs  at  the  data  frequency  yf,  .  The  frequency  component  of 
at  0.  29  y^O  is  20  db  attenuated  by  this  function. 


-1(10) 
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Consequently,  data  analyzed  in  the  manner  outlined  can  be 
used  only  over  a  very  small  frequency  interval  near  the  maximum 
frequency  components.  Because  the  magnitude  of  the  random  phase 
fluctuations  at  frequencies  near  the  interferometer  fringe  frequency  is 
expected  to  be  negligibly  small,  the  output  from  this  data  program  is 
not  considered  to  be  satisfacory,  and  the  technique  was  discarded. 
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Appendix  C>2 

« 

Accuracy  Requirements  for  Aircraft  Positions  Data 


The  phase  difference  between  two  antennas  of  an  inter¬ 

ferometer  of  baseline  (c|)  is  given  (when  parallax  is  neglected)  by 


4m  ^  ^  ‘f'v 

/S 


C-2(i) 


where  ^  is  the  unit  direction  vector  from  the  antennas  to  the  source, 
and  is  the  total  random  phase  noise  in  the  system,  which  inclid  es 
phase  path  fluctuations,  source  vibrations,  thermal  noise,  etc.  If 
were  zero,  then  the  measured  phase  difference  would  be  exactly 
equal  to  ^.(3  )  .  In  the  presence  of  phase  noise,  however, 

we  will  make  an  error  ^  in  the  determination  of  the  position  vector 
BO  that 


•d  i-  S  *  cl  J 

A  ^ 


C-2(2) 


'S 

where  /C  is  the  true  unit  position  vector  and  is  the  "unit  angular 

error"  vector,  i.e.  /\  is  a  vector  normal  to  ,  so  that  f  | 

is  the  angular  error  in  the  determination  of  ,  The  error  term 

in  Equation  C-2(2)  is  ? Ifd ^  ^  and  this  term  must  be  equal  to, 
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or  omaller  than 


./ai^ 


.  Contoquently  we  can  eetimate 


an  upper  bound  for  the  required  eo  that  maybe 

determined  to  an  accuracy  .  Then 


1^1  - 


C-2(3) 


In  our  experiment^  we  wieh  to  determine 


to  an  accuracy 

riment 


of  0.5^  10*^  radian*.  Also  for  our  experiment 

X 

that 


radians 
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Comparison  of  Polynomial  Coefficients  to  be 
Expected  from  Aircraft  Position  Data  with 
Coefficients  of  Fit  to  the  Interferometer  Phase  Data  * 

In  our  discussion  of  results  we  have  implied  that  the  coefficients 
of  fit  B«  representative  of  the  coefficients  ^ t  (a,  and  c  in 

the  polynomial  that  describes  aircraft  motion.  To  verify  this  contention, 
we  calculated  the  linear,  quadratic  and  cubic  coefficients  in  equation  (C-lOc) 
from  aircraft  tracking  data  recorded  during  the  flight  test  on  28  August  1963. 
These  data,  of  course,  are  not  sufficiency  accurate  for  determining 
aircraft  position  to  the  extent  required  for  this  experiment,  but  are 
adequately  precise  to  illustrate  the  magnitudes  of  the  coefficients 
involved  here. 

On  28  August  1963,  at  16^  23*^  40  sec  EDT,  the  initiaUposition 
vector  had  a  magnitude  of  35  N  miles  and  an  inclination  to 

the  north-south  interferometer-baseline  of  110°.  The  aircraft  velocity 
vector  ( V)  had  a  magnitude  of  428  kts  and  an  inclination  to  the  north-south 
interferometer  baseline  cf  191°  (see  diagram).  We  estimate  that  the 
confidence  attached  to  these  position  data  is  as  follows; 

a)  velocity  error  of  ±10% 

b)  aircraft  course  relative  to  baseline  ±7° 

c)  range  error  of  ±7% 

d)  angle  between  baseline  and  position  vector  ±0.  5° 
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F rom  these  data  we  calculate  a  linear  term  (equal  to  O  )* 

144  rad/sec,  a  quadratic  term  (equal  to  C>  )  of  0.17  rad/sec^ 

and  a  cubic  term  (equal  to  ^  )  of  6. 3  x  10*^  rad/sec^.  The 

linear  and  quadratic  terms  (  0  and  in  equation  (C- 13c)) obtained 

from  a  least>squar«s  fit  to  5  seconds  of  analog  phase  data  (from  16^  23"* 

38  sec  to  16^  23*"  4.'  sec)  were  fotmd  to  be  164  rad/sec  and  0.16  rad/sec^, 

respectively.  A  cubic  term  was  not  computed  in  our  machine  program; 

however,  it  can  be  shown  that,  due  to  inadequate  knowledge  of  the  exact 

aircraft  position  the  cubic  term  ^ (L  computed  above  will  be  in 

-4  3  ^ 

error  by  ±2.2  x  10  rad/sec  .  In  Figure  C-3-2  we  have  plotted  the 

2 

spectrum  of  JT along  with  the  spectrum  of  ZJL  • 

*'  A 

The  two  spectra  are  very  similar  to  one  another,  particularly  when  the 
limitations  on-  included.  However,  the  uncertainty  in 

determining  aircraft  velocity,  course,  and  range  is  too  great  to  state 
positively  that  the  spectrum  of  "JT is  primarily  due  to  either 
aircraft  motion  or  to  atmospheric  effects. 


s 


See  Equation  (C»8). 
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Figure  0-^-1  AIRCRAFT  POSITION  DATA  FOR  28  AUGUST  1963 
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Appendix  C-4 

Weather  Summaries  and  Aircraft  Position  Data 

August  1,  1963  (Figure  C-4-1) 

A  low  pressure  trough  passed  Western  New  York  during  the 
morning  hours  inducing  west  to  west-northwest  winds  at  the  surface 
and  aloft  during  the  afternoon.  Showers  and  thunderstorms  dominated 
the  weather  picture  over  Western  New  York,  Southern  Ontario,  and 
Lakes  Erie  and  Ontario  all  afternoon.  These  storms  moved  generally 
eastward  to  east-southwestward  at  about  30  mph. 

At  1300  EST  a  line  of  storms  was  oriented  from  Rochester,  N.Y. 
south- southwest  through  Olean,  N.Y.  with  the  strongest  cells  near 
Olean.  This  line  of  storms  continued  eastward  to  lie  along  a  line 
from  about  Auburn,  N.Y.  to  Wellsville,  N.Y.  by  1545  EST. 

Thunderstorms  covered  most  of  the  area  north  and  northwest 
of  Buffalo  ail  afternoon.  The  southern  edge  of  this  group  of  storms 
extended  from  Port  Colbourne,  Ontario  to  Lyndonville,  N.  Y.  at  1300  EST, 
At  this  time  the  strongest  cells  were  located  near  Toronto,  Ontario. 

The  southern  edge  of  the  area  of  thunderstorms  extended  from  Hilton, 

N.Y.  west  to  Welland,  Ontario  at  1345  EST  and  extended  south  to 
State  Highway  33  by  1445  EST.  By  1545  EST  the  southern  extreme  of 
the  thunderstorms  was  about  where  it  had  been  two  hours  earlier. 
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Figure  0.4-1  AIRCRAFT  POSITION  DATA  FOR  I  AUGUST  1963 
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Throughout  this  two-hour  period  the  heaviest  activity  continued  along 
a  line  south  of  Lake  Ontario  from  Hilton,  N.  Y.  west  almost  to  Niagara 
Falls. 


Scattered  thunderstorms  and  showers  persisted  most  of  the 
afternoon  over  Eastern  Lake  Erie.  Little  or  no  activity  occurred  to 
the  southwest  toward  Jamestown  during  the  afternoon. 

August  22,  1963  {Figure  C-4-2) 

The  same  high  pressure  system  that  dominated  the  weather 
pattern  on  the  21st  persisted  but  weakened.  A  few  sc~>ttered  thunder¬ 
showers  developed  over  Southern  Ontario  prior  to  1545  EST.  At 
1545  EST  a  tine  of  showers  extended  from  Bowmansville,  N.  Y. 
northwestward  to  Shelbourne,  Ontario  and  an  isolated  thundershower 
was  located  about  ten  milss  west  of  Niagara  Falls.  These  storms 
moved  eastward  at  about  15  mph  and  dissipated.  By  1745  EST  there 
were  only  a  few  showers  in  the  vicinity  of  Toronto,  Ontario,  and  a 
weak  line  of  showers  north  of  Kingston,  Ontario  between  Peterborough, 
Ontario  and  Ogdensburg,  N.  Y. 
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Figure 


AIRCRAFT  POSITION  DATA  FOR  22 


AUGUST  1963 
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August  27,  1963  (Figure  C-4-3) 

A  Continental  Polar  high  pressure  system  was  centered  over 
Western  New  York  supressing  the  formation  of  clouds  and  precipitation. 
No  significant  precipitation  was  observed  within  100  miles  of  Buffalo. 

August  28,  1965  (Figure  C-4-4) 

The  eastward  retreat  of  a  Continental  Polar  high  pressure 
system  and  the  approach  of  a  low  from  the  southwest  coupled  to 
produce  wide  spread  showers  over  Western  New  York  State  during  <  e 
afternoon  and  evening.  At  1050  EST,  rain  showers  were  reported 
across  Southern  Ontario  moving  eastward  at  about  30  mph. 

By  1450  EST  a  large  area  of  precipitation  aloft  was  spreading 
across  Western  New  York  State  and  into  Central  New  York.  Dry 
surface  air  prevented  much  of  the  precipitation  from  reaching  the 
ground  until  about  1600  EST.  By  1645  EST  light  showers  generally 
extended  to  the  ground  and  covered  the  area  northwest  of  a  line  from 
Rochester,  N.  Y.  to  Erie,  Pa.  Somewhat  heavier  showers  were  then 
occurring  over  southern  Ontario.  F'recipitation  occurred  at  Buffalo 
between  1232  EST  and  1248  EST,  and  from  1610  EST  until  1910  EST. 
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Fisure  C-M  OIRCMFT  POSITION  ORU  FOR  27  AUGUST  1963 
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Figure  AIRCRAFT  POSITION  DATA  FOR  28  AUGUST  1963 
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APPENDIX  D 

Atmospheric  Propagation  Effects  on  a 
Focussed,  Multiple  Antenna,  Communications  System 


A  complex  of  four  thirty-foot  X-band  antennas  has  been  proposed 
for  the  terminal  equipment  of  a  satellite  communication  system.  The 
four  antennas,  positioned  at  the  corners  of  a  square  10  miles  on  a 
side,  will  be  phased  to  obtain  coherent  illumination  of  the  satellite. 

A  passive  or  semi-active  satellite  at  an  altitude  of  2200  miles  is  envisioned 
and  communication  between  ground  terminals  is  desired  whenever  the 
satellite  is  more  than  six  degrees  above  the  horizon  of  both  terminals. 

The  purpose  of  this  document  is  to  indicate  the  nature  of  the 
propagation  problem  and  show  how  the  proposed  system  will  be  affected 
by  propagation  anomalies. 

Ground  Rules  of  Study 

For  the  purpose  of  this  study  it  is  assumed  that  the  effective  gain 
of  a  60  foot  parabola  is  required  at  X-band  (9000  Mcs).  If  four  thirty-foot 
antennas  are  to  be  used  to  achieve  full  gain,  they  must  be  continuously 
phased  so  that  the  satellite  "rides"  the  maximum  of  the  interference 
pattern  -  the  phase  of  each  of  the  four  transmitters  must  be  such  that 
the  individual  fields  at  the  satellite  ere  m  phase.  We  shall  call  this 
phasing  operation  "focusing"  the  antenna  system. 
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It  should  be  noted  that  if  a  single  thirty«toot  parabola  for 
transmitting  and  receiving  were  used  in  place  of  a  60*foot  antenna, 
the  system  loss  would  be  increased  by  12  db.  If  the  transmitter  power 
is  divided  equally  among  four  30  foot  transmitting  antennas,  and  four 
phased  30  foot  antennas  are  used  for  receptioi*  but  the  pliasing  control 
is  so  poor  that  ah  phases  (between  0  and  2ir  radians)  are  equally 
probable  for  the  transmitting  and  receiving  complexes,  then  not  only 
will  the  system  loss  be  increased  by  12  db  but  the  received  signal 
will  also  exhibit  deep  and  quasi-Rayleigh  fading.  The  "Rayleigh" 
fading  will  degrade  the  resultant  signal  to  the  point  that  it  would  be 
better  to  transmit  foil  power  into  a  single  30  foot  dish,  and  receive  with 
a  single  30  foot  dish.  Overall  phase  control  of  the  transmitting  and 
receiving  antennas  is  clearly  important. 

Propagation  Effects 

The  refractive  index  of  the  earth's  atmosphere  is  a  point  by 
point  random  variable.  Experiments  conducted  by  the  National  Bureau 
of  Standards,  among  others,  indicate  that  even  on  a  fixed  geometrical 
path,  the  phase  of  the  received  signal  is  a  random  variable.  When 
transmissions  from  a  single  antenna  are  received  at  two  other  separated 
antennas  and  ti  e  instantaneous  phase  'ifference  between  the  two  received 
signals  is  recorded,  the  phase  difference  is  found  to  be  a  random  variable. 
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Figure  D-1  is  a  plot  o£  the  power  spectrum  of  the  random  phase 
difference  measured  between  antennas  separated  by  approximately 
one  mile  (NBS  Hawaii  Data).  The  length  of  the  propagation  path 

was  15.  5  iniles  with  the  transmitter  situated  on  a  10,000  ft.  peak  and 
the  receivers  located  at  an  altitude  of  100  feet  along  a  baseline  perpen*** 
dicular  to  the  prop  .gation  path.  The  spectral  density  estimates  shown 
are  based  on  15  minute  sample  lengths  and  the  spread  in  the  ordinate, 
value  ■>  is  indicative  of  the  variability  of  the  data. 


Figure  D-2  is  a  plot  of  the  power  spectrum  of  the  phase  fluctuation 
measured  along  a  single  path  during  the  NBS  Hawaii  program  (Ref.  1). 

If  4^^  is  the  variance  of  the  phase  fluctuations  measured  along  a 
single  path  and  ^  is  the  variance  of  the  difference  of  phase  between 
two  antennas,  then  one  would  expect. 


do-  ^ 


(EM) 


where  ^  is  the  normalized  cross  correlation  coefficient  of  the  phase 
fluctuations  at  the  two  antennas. 


As  the  separation  between  the  antennas  increases,  one  would 
e}cpect  the  variance  of  the  phase  fluctuation  data  to  increase  until  the 
limiting  value  of  ^  was  reached. 
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CYaES  PER  SECOND 

Figure  D-l  HALEAKALA  -  PUUNENE  PATH:  POWER  SPECTRA  OF  Mc/s 

PHASE  DIFFERENCE  VARIATIONS.  BASELINE  LENGTH  4847  FT(A-H) 
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Figure  D-2  HALEAKALA  -  PUUNENE  PATH:  POWER  SPECTRA  OF  Mc/s 

SINGLE-PATH  PHASE  VARIATIONS 
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Figure  D-3  is  a  plot  of  the  variance  of  the  phase  difference 
between  two  antennas,  plotted  as  a  function  of  base  length  (also  NBS 
data  Ref.  1).  There  is  a  general  increase  in  the  variance  as  the 
baseline  length  is  increased  until  at  something  like  1000  feet  the  curve 
levels  off  near  the  single  path  value.  It  should  be  noted  that  the  data 
presented  are  not  simultaneous  (which  accounts  for  the  large  spread  of 
variances)  and  Figure  D<-3  is  presented  only  to  show  the  trend  of  the  data. 

Figure  D-4  is  a  compilation  the  available  single  path  phase 
fluctuation  data  (Thompson,  Janes  &  Kirkpatrick  I960  )  normalized 
for  frequency  and  path  length. 


It  should  be  noted  that  the  median  Hawaii  data  indicate  spectral 
densities  an  order  of  magnitude  greater  than  the  median  of  the  Colorado 
data.  It  should  be  further  noted  that  during  a  period  in  Colorado  when 
a  frontal  rain  storm  occurred  over  the  propagation  path,  an  apparent 
increase  in  spectral  density  above  the  median  density  of  some  20  db 
was  observed  and  these  data  are  indicated  by  the  circled  data  points. 

As  a  model  of  the  troposphere,  assume  that  the  troposphere 
extends  to  30,000  feet.  If  we  take  as  a  reference  tropospheric  path,  a 
path  at  an  elevation  an^  of  six  degrees  extending  from  ground  level  to 
30,000  feet,  our  reference  tropospheric  path  length  will  be  300,000  feet 
or  3,000,000  wavelengths  at  X-band.  The  ordinate  value  of  Figure  D-4 
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BASELINE  LENGTH  IN  FEET 


gure  0-3  HALEAKALA-PUUNENE  PATH:  VARIANCE  OF  PHASE  DIFFERENCE 

VERSUS  BASELINE  LENGTH 
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Figure  D>4  SPECTRUM  OF  APPARENT  PATH  LENGTH  VARIATIONS 
DUE  TO  THOMPSON.  JANES  AND  K I RKPATRI CK 


D-M 


can  then  be  renormalized  to  our  reference  path  by  multiplying  the 
ordinate  by  355  (2v  x  3)^  to  convert  it  to  square  radians/unit  bandwidth. 
This  has  been  done  in  Figure  D-5. 

Figure  D-5  should  now  be  viewed  as  the  long  term  power 
speccrum  of  the  range  fluct'iations  to  be  expected  on  a  single  path  for 
a  fixed  transmitter  and  receiver  location.  We  will  first  derive  the 
phase  difference  spectrum  fur  an  interferometer  with  a  ten  mile  baseline. 
If  one  identifies  the  temporal  fluctuations  of  phase  on  a  single  path  with 
the  drift  of  a  turbulent  refractive  field  across  the  propagation  path, 
it  is  possible  to  use  single  path  spectra  to  deduce  the  phase  difference 
spectra  as  follows: 

A  wind  speed  of  10  feet  per  second  (appropriate  to  Maui)  would 
cause  spatial  irregularities  of  10  feet  to  create  a  one  cycle  per  second 
fluctuation.  A  scale  relating  frequency  to  irregularity  size  is  indicated 
in  Figure  D- 6.  Also  on  Figure  D-6  is  the  estimate  of  spectral  density 
we  shall  use  in  our  calculations.  As  can  be  seen,  it  is  approximately 
25  db  more  pessimistic  than  the  median  Colorado  spectral  density  and 
approximately  5  db  more  pessimistic  then  the  Colorado  rain  data.  It 
is  expected  that  our  estimates  cf  system  performance  will  therefore 
reflect  the  worst  possible  performance.  If  the  system  can  perform 
under  these  circumstances,  it  should  perform  well  most  of  the  time. 
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Figure  [>-5  EXPECTED  POU/ER  SPECTRUM  OF  PHASE  FLUCTUATIONS 
AT  X-BANO  FOR  A  300,000  FT  TROPOSPHERIC  PATH 


I 


The  spectrum  of  Figure  D>6  has  been  further  modified  to  account 
for  the  finite  baseline  length.  It  is  clear  that  irregularities  smaller 
than  10  miles  will  cause  uncorrelated  phase  fluctuations  at  the  two 
antennas.  Irregularities  much  larger  than  10  miles  would  cause 
correlated  phase  fluctuations  at  both  antennas  such  that  the  phase 
difference  between  the  two  antennas  would  be  zero.  It  can  be  shown 
that  the  effect  of  the  finite  baseline  can  be  accounted  for  by  multiplying 
the  single  path  spectrum  with  a  filter  function  such  that  all  spatial 
frequencies  greater  than  twice  the  baseline  length  are  attenuated  with 
an  f^  slope.  The  effect  of  the  finite  baseline  is  indicated  in  Figure  D-6 
by  the  labelled  dotted  line. 


We  expect  that  the  major  contribution  to  the  fluctuation  spectrum 
will  be  dominated  by  the  altitude  region  2000-5000  feet.  A  wind  speed 
of  20  feet  per  second  appears  to  be  representative  of  this  altitude. 

Accordir.gly,  the  spectrum  in  Figure  D-6  has  been  modified  to  i.\ke 
account  of  a  20  fps  wind  speed  by  multiplying  the  abscissa  values  bv  2  and 
dividing  the  ordinate  values  by  2.  The  resulting  spectrum  is  presented  in 
Figure  7  and  represents  the  expected  rate  of  the  spectrum  of  phase  difference 
fluctuations  experienced  by  a  10-mile  baseline  system  operating  in  the  U..  S. 
at  X-band  using  30-foot  antennas.  It  should  be  noted  that  this  spectrum  strictly 
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Fi guru  D-7  SPECTRUM  OF  PHASE  DIFFERENCE  FLUCTUATIONS  EXPECTED  IN  U.  S. 

FOR  10  MILE,  X’BANO  BASELINE  USING  30  FT  DIAMETER  ANTENNAS 
300.000  FT  PATH  AND  MEAN  WIND  OF  20  FPS 


represents  the  situation  when  the  target:  is  fixed.  If  the  target  moves, 
the  line  of  sight  is  rotated  through  a  turbulent  field  and  <xie  might 
expect  the  spectrum  to  have  higher  frequency  components  than  does 
the  spectrum  for  a  fixed  target.  However  in  the  case  of  the  satellite, 
the  linear  velocity  of  the  line  of  sight  evaluated  at  2000-5000  feet  altitude 
is  much  less  than  th^  assumed  mean  wind  speed  of  20  fps  so  that  the 
spectrum  of  Figure  D-7  should  represent  to  a  high  degree  of  accuracy 
the  spectrum  due  to  the  moving  line  of  sight. 


Effects  of  Random  Phase  Fluctuations  on  Power  Density 

Consider  that  the  fields  radiated  from  four  separated  antennas 
are  so  phased  that  the  fields  from  each  antenna  are  in  phase  except  for 
a  random  fluctuation  of  phase  in  each  field  = 

We  shall  assume  that  the  are  independent  chance  variables 
having  a  normal  distribution  with  mean  value  zero  and  standard  deviation 


The  net  field  at  the  focal  ppot  can  be  represented  by 
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The  probability  of  obeerving  a  valtae 


(Dt3) 


and  aince 


Ht. 


are  independent 


Then  the  mean  aquared  value  of  the  field. 


focal  apot  and  in  Appendix  D-1,  it  ia  ahown  that 


,  ia  the  power  in  the 

W  u 


q-  -(-iX  JW 
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This  formula  has  the  expected  variation  with  Q  .As  the 
phase  fluctuations  go  to  zero,  the  power  goes  to  16,  i.e.,  to  4^,  as 
expected.  When  the  phase  fluctuations  become  so  large  that  the  four 
antennas  act  as  four  incoherent  sources,  the  net  power  is  the  sum  of  the 
powers  radiated  (four  in  this  case). 


Consider  a  system  in  which  the  appropriate  phase  corrections 
to  obtain  a  focus  in  the  absence  of  atmospheric  refractive  fluctuations 
were  developed  from  geometry.  Although  it  is  highly  unrealistic  to  do 
so,  we  shall  further  assume  that  the  phase  correction  so  obtained  is 
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perfect  and  we  ihall  then  calculate  the  effect#  of  the  random  nature 
of  the  troposphere  upon  the  power  density  in  the  focal  spot. 

From  Figure  D-7.  .  M| «  10^  .  1. 42  x  10^  per  peth 

Consequently. 

jtT-  ^  ^ 

or  the  system  is  completely  incoherent  so  that  some  sort  of  active  phase 
compensation  is  required. 


Prediction  Problem 

The  previous  section  considered  the  degradation  caused  by 
tropospheric  phase  fluctuations  for  a  stationary  satellite  of  known 
position,  and  indicated  that  phase  fluctuations  would  cause  the  system 
to  fail  u'.less  active  phase  compensation  were  utilized.  In  the  actual  • 
satellite  problem  however,  the  target  is  moving  and  we  are  required 
to  predict,  from  our  measurements,  what  the  required  phase  will  be 
to  focus  at  the  satellites'  new  position. 


One  way  of  accomplishing  the  prediction  problem  is  to  measure 
the  phase  difference  between  antennas  as  a  function  of  time.  If  L 
is  the  travel  time  between  the  satellite  and  the  ground  station  and 
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is  the  phase  difference  between  the  two  antennas  at  a  time  T,  •  the 
phase  difference  at  a  time  ^  ** 


cjxt.U 


(0-^6) 


where  is  the  time  derivative  of  (p  evaluated  at  t  -  ti  . 

The  phase  correction  required  to  lead  the  target  will  be  the  conjugate 
of  Equation (IV6). If  the  atmosphere  were  ideal  and  there  were  no  fluctuations 
of  refractive  index  all  the  variablesti(|>4|  would  be  proportional  to 
the  time  variation  of  geometrical  quantities.  Let  us  denote  the  phase 
factors  appropriate  to  pure  geometrical  factors  with  the  subscript 
zero.  Thus 


(D47) 


where  random  component  of  the  phaae  diffe'cnce  between 

the  two  antennas. 


Combining  and  we  see  the  predicted  correction. 
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(D-(8) 
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Af  will  be  •howii,  the  second  term  on  the  right  hand  side  of  (D-8)  contains 
a  source  of  error  which  must  be  evaluated.  Referring  to  the  basic 
concepts  inherent h(D^,  while  the  rate  of  change  of  the  phase  due  to 
geometrical  considerations  may  lead  to  a  reasonable  prediction,  the 
instantaneous  frequency  (or  rate  of  change  of  phase)  due  to  the  random 
phase  fluctuations  mev  only  lead  to  reasonable  predictions  if  the  lead 
time  is  short  enough.  The  error  in  Equation  (Delias  in  the  second 
term  which  is  a  stochastic  variable.  Define 


and  are  normally  distributed  variables.  The  variance 

of  is 
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In  order  to  evaluate (£>*10)  we  require  the  variance  of  €<*)■  It  can  be 
shown  that  if  4'*)  has  a  power  spectrum  5M 

,  the  power  spectrum  of 

>•  The  variance  of  U)  is  the  area  under 

the  curve  of  the  power  spectrum  and  is 
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The  ordinate  values  of  Figure  D>7  have  been  multiplied  by 


4«^  and  by  the  square  of  the  value  of  the  abscissa  and  the  result  plotted 

as  the  power  spectrum  of  /i)  in  Figure  1X8.  The  area  under  the  curve,  the 
variance  of  ^  is  found  to  be  3.45  (radians /sec)*.  The  standard 

deviation  of  ihe  instantaneous  frequency  fluctuation  is  i.86  radians/eec. 


Now  the  lead  time  required  for  a  satellite  at  2200  miles  and  six 
degrees  above  the  horizon  is  0. 118  seconds.  (This  assumes  that  all 
data  processing  on  the  ground  takes  place  in  a  few  tens  of  microseconds.) 

(DM2) 


On  the  average,  then, we  would  predict  a  change  in  phase,  due  to 
the  random  nature  of  the  atmosphere,  which  has  a  variance  of  .048 
(radians)*.  Now  the  prediction  time  m*.8t  be  the  travel  time  from  the  earth 
to  the  satellite  -  but  the  troposphere  extends  to  only  30,000  feet  and 
the  radio  wave  spends  something  less  than  one>half  millisecond  in  the 
random  troposphere.  We  would  like  to  know  how  much  the  phase  (difference) 
can  change  in  a  half  millisecond. 
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In  Appendix  D>  2  we  show  that  the  mean  squared  phaee  fluctuation 


which  can  occur  in  a  time  i  ia 


^  a)]  ( .  -  p  ^ 


The  mean  square  value  of  the  phase  change  in  one  millisecoml 


is  therefore 


/Oa. 


l~':  S' 


1  <nr‘’i' 3.40  ~  10 


(D-13) 


However,  our  prediction  technique  programs  a  mean  squared  phase 
correction  of  .048  (radians)^. 


The  difference  between  equations  (D-13)  and  (D-1^)  is  an  estimate 
of  the  mean  squared  error  incurred  in  this  prediction  process.  Mean 
squared  phase  error 

(f  Aua  [  -  ora/v  [  A(f  ]  ^  i  ^ 


(D-14) 


The  average  power  in  the  focal  spot  will  then  be 


jEl'-  ^  It 


(D-15) 
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and  this  four  antenna  system  will  act  like  a  coherent  system. 

Satellite  Veloc ity  Correction 

The  previous  analyses  are  approximate  since  they  assume  that  the 
velocity  of  the  satellite  can  be  neglected  in  determining  the  phase  differ¬ 
ence  between  antennas  at  any  instant  of  time. 


T 


The  above  sketch  exaggerates  the  real  situation  but  does  clearly 
illustrat®  the  physical  situation.  The  satellite  at  position  1  is  further 
from  antenna  A  than  it  is  from  antenna  B.  By  the  time  the  signal  emitted 
at  position  1  reaches  antenna  A,  a  signal  emitted  from  position  2, 
travelling  over  the  shorter  path  r2  will  reach  antenna  B.  Satellite 
position  2  is  displaced  a  distance  /V'l  (where  AT  is  the  satellite 
velocity)  from  position  1.  The  time  of  propagation  between  position  1 
and  antenna  A  is  'i 


I 
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The  phase  difference  between  antennas  will  therefore  be 

©  -  UjT 

We  solve  for  ^  as  follows: 

(2 ■* [’'•  + 1  c t *  X .'-t “^4^  =  C^o-*d*'rCl  i'j. 


Let  7  c?  ^ 

C 


and  rearranging  terms  we  get 


/JCi'.t;*')]  -  Xnx^n^ik-O 
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We  use  the  binomial  expansion  of  keeping  all  terms  with 

/x*'"  / 

coefficients  less  than  or  —  to  obtain: 


I. 


The  phase  difference  between  antennas  is  ^ 


u,T>  ‘iii- 


(D-16) 


It  should  be  noted  that  when  the  velocity  of  the  satellite  is 
neglected  (  (^  '^  0*) 


'!Is\ 

T',  L  r'v 


A 


and  the  usual  interferometer  formula  obtains  (far  field  approximation). 

In  order  to  estimate  the  magnitude  of  the  correction  factor  — ZT" 

fXm 

note  that  the  velocity  term  is  maximized  for  large  Xo  implying 
low  elevation  angles.  At  low  elevation  angles  for  a  Z200  mile  satellite 


V,  cl  ^  ^/i  ~  co^ 
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where  E  is  the  elevation  angle  referred  to  the  center  of  the  array. 
Accordingly  we  write  (P*16)  as 


(D-t7) 


where  /V^  is  the  radial  component  of  the  satellite  velocity. 

An  estimate  of  the  magnitude  of  the  velocity  dependent  term 
can  be  made  by  assuming  that  after  the  satellite  emits  its  pulse  at 
position  Z,  we  freeze  the  satellite  in  space  (a  highly  theoretical  proposition). 
We  then  ask,  if  we  transmitted  the  conjugate  of  the  phase,  whether 
we  would  focus  on  the  satellite. 


At  low  elevation  angles 


A 


opt 


uvio^ 


!o 


/V 


*7  ^ 

$=■  ■■ 

3xi6‘' 


The  correction  term  is  of  the  order  of  84  radians  and  since, 
in  general,  the  correction  term  will  not  be  an  integral  multiple  of  Zn 
radians,  it  must  be  accounted  for  if  deep  fading  signals  are  to  be  avoided. 
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Since  the  required  correction  is  proportional  to  d,  the  distance 
between  antennas,  decreasing  the  separation  between  antennas  as 
well  as  increasing  the  minimunn  elevation  angle  will  lessen  the  problem. 
For  example,  an  antenna  separation  of  500  feet  and  elevation  angle 
of  45°  will  result  in  a  required  correction  of  approximately  1/2 
radian  which  may  be  acceptable.  We  conclude,  however,  that  an  active 
phase  compensation  system,  with  an  added  satellite  velocity  correction, 
will  be  required  to  obtain  a  fully  coherent,  efficient  communication 
system  using  widely  spaced  anternas. 

Summary  and  Recommendations 

In  Appendix  D-3  we  examine-  the  accuracy  of  a  linear  prediction 
technique  and  find  that  errors  due  to  linear  prediction  are  negligibly 
small.  In  general  then  we  have  shown  that  using  "typical"  propagation 
data,  the  X-band  communication  system  under  consideration  can  be 
used  down  to  elevation  angles  as  low  as  6°  if  an  appropriate  correction 
for  the  satellite  velocity  is  added  to  the  active  phase  compensation 
mode  of  operation. 

It  should  be  recognized  that  this  fully  coherent  system  at  low 
elevation  angles  has  only  a  12  db  advantage  over  a  far  simpler  system 
employing  single  thirty  foot  dishes  at  each  terminal.  The  cost  of  the 
fully  coherent  system  is  not  known  but  ultimately  the  question  must  be 
answered  in  terms  of  just  how  much  12  db  is  worth. 
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Up  to  this  point,  we  have  not  considered  the  question  of  system 
bandwidth.  It  is  clear,  because  of  the  separation  between  antennas,  that 
in  addition  to  the  radio  frequency  phasing  required  to  focus  the  antenna 
system,  envelope  delay  must  be  inserted  in  the  system.  Without 

envelope  delay,  the  maximum  allowable  bit  rate  would  be  something 

3 

like  3.57  x  10  .  The  maximum  multipath  delay  for  a  14  mile  baseline 
is  140  microseconds.  If  the  data  pulse  from  each  transmitter  is  twice 
that  length,  the  transmissions  from  each  antenna  will  overlap  and 
coherent  operation  can  be  obtained.  While  the  bit  rate  is  greater  than 
that  normally  available  via  conventional  HF  circuits  it  is  still  a  far 
cry  from  the  minimum  desirable  information  bandwidth  of  two  megacycles. 
Envelope  delay  is  a  real  requirement  then.  If  bandwidths  of  greater 
than  two  megacycles  are  to  be  achieved,  the  range  from  each  transmitter 
to  the  satellite  must  be  known  to  better  than  one-half  microsecond. 

The  ability  to  measure  the  range  of  each  antenna  to  the  satellite  and 
the  adjust  envelope  delay  will  set  the  limits  on  the  maximum  bandwidth 
of  the  system. 
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APPENDIX  D-1 


Calculation  of  Power  in  Focal  Spot 


Let  the  fields  radiated  by  each  of  four  antennas  have  unit 
amplitude.  Let  each  field  have  a  phase  factor  which  is  the  sum  of  two 
terms: 


1)  a  term  to  insure  perfect  focus  in  the  absence  of  propagation 
phase  fluctuations: 

2)  a  term  which  represents  the  atmospheric  propagation 

fluctuations.  ■ 


The  second  term  is  a  random  variable.  The  available  experimental 
data  show  that  the  are  normally  distributed  and,  if  the  ancennas 

are  separated  by  more  than  100  meters,  that  phase  fluctuations  are 
uncorrelated.  Accordingly, 


^  A/} 


/yti/ub^  u  u 


I 


4 

T 
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DeftAing  the  reeultiiig  field  megnitode  as  R,  we  write 


4  ^ 

I  ^  X.  ^  ^  j 


Let  US  look  at  cosine  terms 


Uo(f. 


because  of  the  independence  of  the  phases 


,r 


r^  (cw -(tws)  ^ 
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The  sine  term  integrals  are  of  form 


r  r 


because  of  the  indep  ndence  of  the  probability  diatribotioos  and 


Therefore 


-  4  4  tl  e 


We  can  readily  generalize  our  result  to  the  power  in  the  focal  spot 


resulting  from  the  phasing  of  n  identical  antennas  and  transmitters 


/y\  h  /h  (ai-»)  e 
« 


For  n  sufficiently  large  and  small 


/y^  €L 
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APPENDIX  D-2 


Rate  of  Phase  FI  actuation 

If  ^)is  a  chance  variable  from  a  Gaassiaa  process  Mich  tiiat 
the  probability  of  obtaining  a  value  between  and  i.  p(*)df 

then 

/  L  f  7 
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The  probability  that  can  change  by  an  aunount  A<P  i 
a  time  'c  i* 

'  -  4- Ml 

l  4f' 


m 


p(M)- 


where  ^  is  the  variance  of 

^fr)is  the  normalized  autocorrelation  function  of  (j> 

It  is  clear  that  A<})  is  normally  distributed  with  variance 

2  f  «>]  --  [=  ['  -  f 


D-2(2) 


D-2(3) 
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APPENIXX  D-3 


Linear  Prediction  Theory 

In  this  section  we  shall  estimate  the  error  accruing  from 
simple  linear  prediction.  The  basic  problem  arises  from  the  fact 
that  linear  predicti'^n  is  a  first  order  approximation  to  a  power  series 
expansion.  A  function  analytic  in  a  region  may  be  expanded  in 
a  converging  power  series 

Linear  prediction  uses  only  the  first  derivative  of  the  function  and 
assumes  that  the  higher  order  terms  are  negligble. 

Consider  the  geometry  sketched  in  the  accompanying  figure, 
is  the  radius  of  themrth.  r.,  represents  the  radius  of  the 
satellite  orbit.  ^  the  phase  difference  between  two  antennas 

separated  a  distance  d  (using  far>field  approximation) 
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where  =; 


Defining 

£  =  v/>  a^tV  [vfU')+'i''<*^''^]^  vJi'w>T\  I'fOl,*-- 
„  -H.  >-  ^ 

(I,  'U) ,  ^  ^  ^ 

cly  ^9- 
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=  to  cwV  ^ 


V'm  -'  iw'^  ^  -  3  st-x  ‘^li  -  c^pI^  f^j 

L  df  d9^  y^J 


t-kcv}9 

l4.lC'-lkUi6- 


dQ'- 


[  i?\  0  k  j_ 

Q  4-  k"^  -T-  ^ 


We  shall  evaluate  the  error  incurred  in  linear  prediction  theory  by  neglecting 
the  higher  order  terms  of  the  power  series  at  Zenith  angles  78,  and  84  degrees. 


X-  0^ 
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After  some  algebra 


=  -  1.56  X  10"*^ 
'T  ijf'Li')  =  4.38  X  10"^ 

-/«»  ^  /4 

£y=  -7.  f^V/0  f  2.l6if/0  f.  .. 

T  ^  7.72.xi6“^ 

f  imx'o 

Cy 

'Z  uiVo  -  3.4yio"^ 


It  would  appear  therefore  that  linear  prediction  theory  is 
more  than  adequate. 
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APPENDIX  D-4 


Far  Field  Approximaticm 


With  reference  to  Figure  D-4-1  on  page  D-39: 
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exact  phage  difference  between  the  two  antennas  = 

I  f 


where 


^  L. 
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Error  in  caused  by  neglect  of  higher  order  terms  is 
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We  must  compare 


f ,  =  i  ( ’t  f  ^ 
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With 
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